ON THE SPECTRAL PROPERTIES OF L± IN THREE DIMENSIONS 

OVIDIU COSTIN, MIN HUANG, AND WILHELM SCHLAG 

Abstract. This paper is part of the radial asymptotic stability analysis of the ground 
state soliton for either the cubic nonlinear Schrodinger or Klein-Gordon equations in 
three dimensions. We demonstrate by a rigorous method that the linearized scalar op- 
^^ ' erators which arise in this setting, traditionally denoted by L±, satisfy the gap property, 

^^ , at least over the radial functions. This means that the interval (0, 1] does not contain 

any eigenvalues of L± and that the threshold 1 is neither an eigenvalue nor a resonance. 
The gap property is required in order to prove scattering to the ground states for so- 
lutions starting on the center-stable manifold associated with these states. This paper 
therefore provides the final installment in the proof of this scattering property for the 
cubic Klein-Gordon and Schrodinger equations in the radial case, see the recent theory 
Q, ■ of Nakanishi and the third author, as well as the earlier work of the third author and 

,^^ ' Beceanu on NLS. The method developed here is quite general, and applicable to other 

. , spectral problems which arise in the theory of nonlinear equations. 



> 



1. Introduction 
1.1. The nonlinear context. Consider the nonlinear Schrodinger equation 



^; iat^-A^ = ±|V^rV it,x)GRl^/ (1) 

f— N . with powers 1 < p < 2* — 1 where 2* = -^^^ in dimensions d ^ 3 and 2* = oo in dimensions 

, I I d = 1, 2. Assuming that ^(i, x) is a smooth solution of sufficient spatial decay, one verifies 

Q ■ by differentiating under the integral sign that mass and energy are conserved: 

Mm)] := ^wmwl = Mmo)] 

5t , The range p < 2* — 1 is referred to as energy subcritical regime due to the fact that in 

the conserved energy the nonlinear term ||'(/'||p+i is controlled by the i/^-norm of ip{t) via 
Sobolev embedding. 

The choice of sign in front of the nonlinearity of ([1]) is crucial: the — sign (known as the 
defocusing nonlinearity) leads to a positive definite conserved energy and one has global 
existence and scattering to a free wave for any data in H^{W^), see [8], [28], and |29j for an 
account of these classical results. Scattering here refers to the property that there exists 
^0 £ H^ so that with the associated free solution ipo{t) := e~ ipo 

mt)-Mt)\\m^O t^oo 

On the other hand, the focusing sign -|-|^|p~^^ on the nonlinearity renders the energy 
indefinite and finite-time blowup may occur, for example for all data of negative energy 
and finite variance, see Glassey [15]. Blowup here refers to the property that ||''/'(*)IIhi ~^ 
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00 as t — 7- T— < oo. In addition, the focusing nonlinearity admits special stationary wave 
solutions of the form e"**" 4>{x) with a 7^ 0, where 

-aV + A<A=l<Ar"V (2) 

Existence of nontrivial decaying solutions to this equation has been known for a long 
time, see for example Strauss [27] and Berestycki, Lions [5]. On the one-dimensional line, 
there are exactly two nonzero decaying solutions which are given by 

Q(x) = ±acosh"^(/3x), a = (^^)^, /3 = ^^ 

The existence and uniqueness can be read off from the phase-portrait in the ((/>, (/)')-plane. 
In higher dimensions no explicit formulas exist and one obtains existence via variational 
arguments. Moreover, uniqueness in the strong sense as in one dimension fails, as it is 
know that there are infinitely many solutions [S]. However, there is exactly one positive, 
radial solution called the ground state. In fact, any positive decaying solution of Q is 
necessarily radial about some point (by Gidas, Ni, Nirenberg [13]) as well as exponentially 
decaying. This unique solution is called ground state and it is the one we consider in this 
paper. 

The orbital stability analysis of this ground state standing wave was settled many years 
ago and depends on the power of the nonlinearity, see [5], [9], [16], [30], |3T]: iov p < p2 := 

1 + 1 (the latter being called L^-critical power) the ground state is stable, whereas for 
P2 =^ P < 2* — 1 the ground state is unstable in the orbital sense. In fact, the instability 
is very strong: arbitrarily small perturbations of initial data Q with respect to the H^- 
topology can lead to finite-time blowup, see [3], [H]) [IH] (and for Klein-Gordon [20]). The 
transition at the power p2 can be seen at the linearized level. More precisely, linearizing 
about the standing wave with a = 1 (which we may assume by scaling) and splitting into 

real and imaginary parts yields the matrix operator Ti := j ~ where 

— J-j^ 

L_ = -A + l-QP-i, L+ = -A + l-pQP-i 

One then finds that for p ^ P2 the spectrum of Ti lies entirely on the real axis, whereas for 
p > P2 there exists a pair of imaginary simple eigenvalues, signifying exponential linear 
instability. 

The more difficult asymptotic stability problem was considered in Buslaev, Perelman [6], 
[7], Soffer, Weinstein [23], Cuccagna [11], but it would take us too far to review the liter- 
ature on this topic. More relevant for our purposes is the conditional asymptotic stability 
problem which refers to the following question: in the orbitally unstable regime, does the 
ground state remain asymptotically stable in forward time under a finite co-dimension 
condition on the perturbation? In fact, due to the structure of the spectrum in that case 
one might expect that a co-dimension 1 condition should suffice in order to stabilize the 
ground state. This is indeed the case, as shown in the orbital stability sense by Bates and 
Jones m for the nonlinear Klein-Gordon equation and for the NLS equation in |13j (fol- 
lowing Bates and Jones) . [1] implemented the Hadamard (or graph transform) method in 
the infinite-dimensional setting given by nonlinear dispersive Hamiltonian PDEs such as 
NLS and Klein-Gordon. The graph transform together with the Lyapunov-Perron fixed 
point approach constitute the only two known methods available for the construction of 
invariant manifolds, and they were both intensely developed in finite dimensions (in other 
words, for ODEs). See the introduction of [191 and the references cited there. 

The asymptotic stability question for the cubic NLS in three dimensions was studied 
in |22] and |3], where the existence of a center-stable manifold near the ground state was 
established on which the solutions remain asymptotically stable and scatters to the ground 
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state. See |17j for the one-dimensional case, and [19], [25j for the Klein-Gordon equation. 
Finally, in the monograph by Nakanishi and the third author [19] (see the references there 
for the original papers) it was shown that these center-stable manifolds divide a small 
ball around the ground state into two halves which respectively give rise to blowup in 
finite positive time on the one hand, and global existence in forward time and scattering 
to zero, on the other hand. 

The most delicate part of the conditional asymptotic stability analysis turns out to be 
the scattering property of solutions starting on the center-stable manifold. This refers to 
the fact that solutions starting on the manifold decompose into a ground state standing 
wave (with slightly different parameters - this is the phenomenon of modulation) plus a 
free wave plus a term which is o(l) in the energy space as t — )• +oo. By some dispersive 
PDE machinery this is equivalent to the property that the perturbation of the modulated 
standing wave satisfies global dispersive estimates, such as pointwise decay (as in [22]) or 
Strichartz estimates as in [3]. 

The reason that such global dispersive estimates can be considered delicate lies with 
the fact that they appear to require detailed knowledge of the entire spectrum of the lin- 
earized operator, including the behavior of the resolvent at the thresholds of the essential 
spectrum. In [22], [3] one therefore needed to assume the gap property of L± for the cubic 
power nonlinearity (i.e., p = 3) in M'^. As already mentioned before, this refers to the fact 
that L± have no eigenvalues in (0, 1] and that 1 is not a threshold resonance. And finally, 
the gap property which we verify in this paper implies via the Lyapunov-Perron method 
that solutions on the center-stable manifold scatter to the ground state, see |19j . 

Demanet and the third author |12] implemented the Birman-Schwinger method numer- 
ically and showed that this assumption is indeed correct (in the general nonradial setting) . 
Moreover, they found that the gap property is even more delicate than expected: it fails 
if the power p on the nonlinearity is lowered slightly below p = 3. This is somewhat 
surprising, as Krieger and the third author |17j . based on Perelman |21] . had shown by 
analytical arguments that the gap property holds in the entire L^-supercritical regime in 
one dimension. This was facilitated by the explicit form of the ground state in dimen- 
sion 1 and one finds, moreover, that L± retain the gap property for all powers down to 
the completely integrable cubic NLS, where a threshold resonance appears. 

While [12] appears to be accurate on all empirical accounts, the numerical method 
implemented there is not a proof since it seems very difficult — if not impossible — to give 
rigorous error bounds for all numerical approximations and calculations required by the 
Birman-Schwinger method. For example, an approximate soliton is computed numerically, 
but without any rigorous bounds on the error introduced by this approximation. 

The purpose of this paper is to offer a completely rigorous, albeit rather computational, 
proof which confirms the gap property of L± in the cubic radial case in M^ . 



1.2. The main equations. Over the radial functions, the equation for the ground state 
reduces to 

-y"{r)--y'{r)+y{r)-y\r)=Q (3) 

r 

By Coffman's theorem [lO] there is a unique, positive decaying solution of Q which is 
smooth on [0, oo). It is denoted by Q and called the ground state. The eigenvalue problem 
over the radial subspace now becomes 

L+u = Au; where L+ = — — ^ ; — hi — 3(5 (4) 

dr'^ r dr 
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and 

L^y = Xy; where L_ = --—r — + 1 - Q^ (5) 

ar"^ r dr 

We show that in the gap [0, 1] the operator L+ has no eigenvalues or resonance, and the 

same is true for L_ on (0, 1]. 

1.3. Technical approach. Eq. ([3]) is likely nonintegrable, and no useful closed form 
representation for Q is known. However, in order to resolve the aforementioned gap 
problem, an exact expression of Q or even the exact values of Q{r) are clearly not required: 
a sufficiently accurate approximation, which we denote by Q, suffices. Unfortunately, 
we have found that the required accuracy for \\Q/Q — l||oo in our problem is on the 
order of 10~^. This is a reflection of the phenomenon seen in [12], namely that the gap 
property is only barely correct. More mathematically, it must mean that L± have complex 
resonances very close to the real axis that become eigenvalues in the gap once the power 
is lowered slightly below p = 3. We remark that this phenomenon may also account for 
the failure of "softer" approaches based on bounds on the number of eigenvalues etc. 

We proceed as follows. We find a suitable approximation Q in the form of a piecewise 
explicit function; for r ^ 5/2 it is given by 

Q{r) = y^{r; /3) = r~Vie~'' + I32g{r); g{r) := r-^(2e''Ei(-4r) - e-''Ei(-2r)) (6) 

for specific /3i,/32, see §2.11 On [0,5/2), the reciprocal 7^ is a piecewise polynomial. The 
coefficients of the polynomials arising in this construction are listed in the first two lines 
of the table on page [3H and we refer the reader to Section [2] for more details of the 
construction, see especially Definition 12.31 

The representation of Q described above is found in the following manner. For r ^ 
5/2 we iterate the Volterra equation once to obtain 7/3. On [0,5/2), we construct by 
Taylor series a solution which is well-behaved at zero and then we extend it by matched 
Taylor series up to 5/2. We take the value <5(0) as a parameter and determine it so 
that Q{b/2~) = 7/3(5/2"''). We then optimize the polynomial representation by sampling 
points from the collection of reciprocals of the aforementioned Taylor series and using 
least squares fitting with three polynomialiJ. We then rationalize the coefficients of the 
polynomials by suitably accurate truncated continued fractions. This is the procedure by 
which the rational numbers listed in the tables at the end of the paper are obtained. 

The next step is to show that Q is close io Q. At this stage the problem is already, in 
some sense, perturbative: 5 = Q — Q \s very small. We thus proceed in a natural way, by 
solving a contractive equation for 5. A slight hurdle arises at this point since the Green's 
function G(r, r') in the integral equation for 5 is not explicit either: G solves a linear 
second order ODE with nontrivial coefficients (combinations of exponential integrals and 
polynomials). We overcome this by finding a nearby equation with explicit solutions, and 
contract out the difference between the two equations. 

Estimating the remainder as a result of replacing Q by Q fi^ © reduces to bounding 
rational functions with rational coefficients. This is done rigorously, as the degrees of 
the polynomials in the denominator and numerator are manageable. There are many 
ways to estimate polynomials. Perhaps the most straightforward one is to place absolute 
values on all coefficients, or at least on all coefficients of powers higher than three, say: 
a polynomial with positive coefficients is easy to bound by monotonicity; it is largest 
at the largest argument. However, inspection of the tables at the end reveals that this 



(i) This is essentially the discrete version of L^([a,6]) orthogonal projection using Legendre polyno- 
mials. Projecting on Chebyshev polynomials would provide an even more economical representation, but 
we prefer the simplicity of least squares fitting; (ii) we look at the reciprocals since they are smoother in 
that the singularities in C are farther away from the real axis and result in more efhcient representations. 



ON THE SPECTRAL PROPERTIES OF L± IN THREE DIMENSIONS 5 

method cannot be applied here, as the coefficients of higher powers are not within the 
smah interval we need (which is < 10~^). In order to overcome this problem, we re-expand 
the polynomials at a number of intermediate points selected so that the coefficients of 
the monomials of degree exceeding 3 are small enough to be discarded modulo small 
errors. Polynomials of degree 3 of course have explicit extrema. The re-expansion refers 
to nothing other than passing from P{r) G Q[r] to the new polynomial in Q[z] defined 
as P{rQ + az), where we always keep \z\ ^ 1. The values of r^ and a are always stated 
explicitly in the text (we use the notation of "partitions" for this purpose). See Note [231 
for more on this issue. 

After having obtained Q up to explicitly controlled errors, we then analyze the spectra 
of L± . This is done essentially in the same way, by finding an accurate Jost quasi-solution 
for r ^ 5/2, and a well-behaved one on [0, 5/2] (which simply refers to the requirement 
that the solution remains bounded with a horizontal tangent at the origin). The quasi- 
solutions are explicit combinations of exponential integrals and polynomials. The way 
the quasi-solutions are obtained also mirrors the soliton approach: iterating the Volterra 
equation for large r and using orthogonally projected Taylor series in the complement. 
This time around, we need the solution ui only on [0, 5/2] and U2 on [5/2, oo]. However, 
since the quasi-solutions depend on the spectral parameter A as well, the calculations are 
more involved. We then check that inf^eio,!] |W^(A)| > (for L+) (inf^eio,!] |W(A)/A| > 
for L~ , resp.) where W is the Wronskian of ui,U2 at 5/2, and this concludes the prooo 

We emphasize that all coefficients involved are in Q, the calculations are exact, and 
the proof, tedious at places - for instance in having to repeatedly re-expand polynomials 
at numerous intermediate points - is fully rigorous. 

In addition, the integral operators upon which the contractive mappings are based have 
small norm (see, e.g., ()14p ). allowing for the calculation of the solutions rapidly and, in 
principle, with arbitrary accuracy. Therefore, this approach is useful numerically as well, 
to obtain rapidly convergent approximants. 

While our approach can in principle be carried out by hand, it is of course unrealistic 
to attempt this in praxis as the calculations in their current form are too numerous as 
well as too long. While we organized Section [2] in such a way that the calculations can 
be in practice done by hand, in later sections we preferred to use the computer algebra 
packages Maple and Mathematica to perform basic operations (such as multiplications of 
polynomials with rational coefficients and solving quadratic equations). The later sections 
involve longer calculations, but of the simple type mentioned above. The exponential in- 
tegral, Ei, is the only transcendental function needed; we estimate it using the asymptotic 
inequalities it obeys and/or by integrating inequalities satisfied by its derivative, an el- 
ementary function. Once more, there are no numerical calculations involved (such as 
numerical integration or numerical location of zeroes) , and with substantially more opti- 
mization effort, it is likely that every step could have been done by hand; we felt there is 
little to gain from this, as human error has a considerably higher chance to occur in such 
a setting. 

We wish to emphasize again that all calculations that were carried out by either Maple 
or Mathematica are completely error free as they only involve finitely many algebraic 
operations in the polynomial ring Q[r]. 

Let us also emphasize that the concrete implementation of the method as it appears 
below is by no means the only possible one. As mentioned in Footnote [T] on the previous 
page, one may substantially reduce the amount of computations required (as well as the 
length of the tables in the appendix) by relying on Chebyschev polynomials instead of least 



A solution bounded near must be a multiple of ui since solutions linearly independent of Ui are 
unbounded at zero; similarly any solution bounded as r — >■ oo is a multiple of U2 . 
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squares fitting in order to carry out the aforementioned projections of the matched Taylor 
series. We intend to present this simpler implementation, together with the nonradial gap 
property in a future publication. 

Finally, we would like to mention that the approach developed in this paper is by no 
means restricted to the specific problem that we study here. For example, in a forthcoming 
publication a similar approach will be used to settle a long-standing spectral question 
arising in completely integrable equations (pertaining to a Painleve equation). 

2. The approximate soliton 

In this section we find an approximation of Q by means of simple functions. Let 

n := 3/10, ra := 17/25, rg := 9/10, r^ := 1, rg := 3/2, rg := 12/5, ry := 5/2 

These increasing numbers define the partition points relative to which we will define the 
piecewise approximations. We denote the characteristic function x([rj,rj)) by Xijd / 7 
and Xi7 = Xi[ri,r7]); similarly x([0,rj)) =: Xoj and x{[rj, oo)) =: Xjoo- We also denote 
Xij = X{[ri,rj]). 

2.1. Solving ([3]) from oo. The following lemma describes all possible solutions of ([3]) 
which decay at oo (at least within a certain range of parameters chosen to suit our needs, 
and up to some error). In what follows, we shall repeatedly encounter the exponential 
integral 



/X gU 
— du; x£ 
-oo u 



where PP denotes the Cauchy principal value of the integral. Define the nonlinear oper- 
ator 

/•oo 

i^if))ir) = / sinh(r - s)s-^f{s) ds (7) 

With gr as in ([6]) one checks that 

rg{r)=M{e-^) (8) 

The exponential integral admits the following asymptotic expansions. 
Lemma 2.1. (i) For each positive integer N one has 



e"^ 



,(-1)'^ . ^.. ^^ .„-.^.,(-l)' 



E^'^<-Ei(-x)<«"^E^!^ Vx>0 (9) 

A:=0 fc=0 

and 

, , e-^'' / 1 3 21 45 \ 15e^ Z'^, ., ,,,,,-,, 



ds 



e-3^ / 1 3 21 45 465 \ 

> ^3~ \^~8 "^ 16^ ~ 64^ ^ 64^ ~ 256^/ 

In particular > g{r) > — e^^'"/8r^ for r ^ rj. 

(a) The function h{r) = —re^g{r) is positive and decreasing. 

(Hi) Define the norm \\ip\\ = sup^j.^^ |^(r)e^| on continuous functions on [r7,oo). Then 

we have 

\\M{e-')\\ ^ 1/8900 (11) 
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Proof, (i) Both ([9]) and (fTO|) follow by means of repeated integrations by parts, (ii) h is 
manifestly positive, while 



h (r) 



.2r 



e-^'s-'^ds<0 



(iii) By (ii), sup^j,^^ |e^AA(e '')| is reached at r = rj and (fTT]) is now immediate from 
(HOD. " D 



Henceforth we assume < /3 ^ 3, which is sufficient for our purposes. Also, recall the 
definition of y^ , see ([6]) : 

2/3(r;/3)=r-i/3e-^ + /3VO 
By the lemma, this defines a positive function for all < /3 ^ 3 and r ^ rj. 

Lemma 2.2. There exists a unique positive solution y{r;(3) to ([3]) with the property that 
y{r; /3) ~ (3r~^e~^ as r ^ oo. It satisfies 



y{r;P) 



2/3(r;/3) 



1 



< 3.2 • 10^^ 



\/ r '^ rj 



(12) 



uniformly in < f3 ^ 3. 



Proof. Setting z{r) = ry[r) (and suppressing the /3-dependence for notational conve- 
nience) yields the ODE 

-z"{r) + z{r) = r^'^z^{r) (13) 

Let ZQ{r) = f3e~^. For solutions z with z{r)e'^ bounded for large r, (J13p can be written as 

z{r) = zo{r)+M{z){r) =: M{z){r) r>0 (14) 

Take the ball B = {h : \\h\\ ^ a/3} and choose a so that AiB C B; the latter condition 
gives 



2^,3 



1 + /3^a 



asCO 



which is satisfied for /3 £ [0, 3] if a = 1 + 1/985, for example. Using (jlip again, expanding 
out (z + 5y — z^ = (5(3z^ + 3z6 + 5^) and estimating each term in the last parenthesis by 
its largest norm in B (such as \\6\\ ^ 2a/3) we get 



\\Miz + 5) - Miz)\\ < 



13a^/3 



2«2| 



8900 



< 



76 



(15) 



Thus A^ is contractive in B and (J14p has a unique solution z^ there. First, since A/'(e~^) = 
©(e"^**) for large r, we have y{r; 13) ~ /3r~^e~'', as claimed. It remains to show (|12p . Note 
that ry'^{r;l3) = Ai{zQ){r). Since 5 := Zs — zq = M{zs) we have 

\\5\\^\\M{zM^P'o?\\M{e~^)\\<^ 
Using this estimate to improve on (llSp we conclude that 



2/(r;/3) 



y3{r;(3) 



Zs-M{zo) 



M{zo) 



^ 1 fJ^a^Mie-nW ... .„_6 
^ 327 /3-/33||AA(e-)|| <'•'■'' 



as claimed. 



(16) 

n 
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2.2. Solving ([3|) on [0,oo) up to a small error. In this section, we study an approx- 
imate solution of ([3]). For the heuristics behind this construction, we refer the reader to 
the introduction. In particular, we specialize the value of (3 in y3(r;/3) in Lemma 12.21 so 
as to most closely approximate the ground state Q{r). 

Definition 2.3. The approximate soliton Q is defined as follows. First, set 

Pi{r) :=gi(r) 

P2{r) := q2{r) + Pi{rs) - g2(?^3) + (Pi(^3) - q'2{r3)){r - r^) ^-^^^ 

r 
^^^"^^ ■" Ae-^ + Brg{r) 

where gi, q2 are the explicit polynomials qi{r) = X^,=q a^ r^ with coefficients as in TableO 
and A^ B are chosen so that P2,P3 match up in a C^ fashion at r = ry; p^ is very close to 
the solution in Lemma 12.21 see ()20p below. Finally, set 

Q ■■= Xos/Pl + X37/P2 + X7oo/P3 (18) 

We remark that with /(r) := ^-^ the coefficients A,B in (J17p are 

^ ^ P2{r7)g' {r-j) +P2{r7)g{r7) 

Pl{r,)ifir-r)g'{rr)-nrr)girj)) 
^. P2{r7)nrr)+p',{n)f{rj) ^ ^ 

P2ir7){f{r7)g'{r7) - f'{r7)g{r7)) 

The specific form of ^3 of course originates with the exact solution from (|12p . From ()19p 
one verifies that 

^<A<^, \B-A'\<33-10^'; B<20 (20) 



We now need to show that Q from (llSp is close to the actual unique ground state Q. We 



begin by checking that Q satisfies ([3]) up to a small error. Below we denote by C^ the 



space of piecewise C^ functions. 

Lemma 2.4. As defined above, Q satisfies the following properties: 
(i) It is decreasing for r G [0,r7], and < Q{r) < 22/5. 
(ii) It belongs to C^{[0,oo)) ilC^ {[0, 00)), andQ'{0) =0. 
(iii) It satisfies the bounds 

Q{r) < 5(1 + r)e-^^x^2{r) + ^e-^''^X27{r) Vr > (21) 

and 

187 e"'' ~, , 350 e"'' , , , 

— — <Qr <— — , forr^rj 22 

69 r 129 r 

(iv) In the complement of the three-point set {r^^r^^r^} the error 

R{r) := -Q"{r) - -Q'{r) + Q(r) - Q{rf 
r 

satisfies the bound 

— 3r 

\R{r)\ < pi(ll/10 - r)xo4ir) + P2(13/5 - r)x47(0 + |^>:7oo(0 (^3) 

6 



where pi := 15-10 and p2 := 25 • 10 
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Note 2.5. (i) To estimate a higher order polynomial P on an interval, we partition the 
interval and Taylor-re-expand P in each subinterval. We define the partitions so that P 
equals the first four terms plus a small error. There are of course other ways to estimate 
polynomials, but this method leads to straightforward calculations. We represent the 
partition by a vector tt, whose components ttj are precisely the partition points. Unless 
otherwise specified, in each interval [7rj,7rj+i) we shall write 

ii{z) = -(1 - z)TTi + -(1 + z)7r,+i, -1 ^ z ^ 1 

and re-expand P{ii{z)) around z = 0. 

(ii) We bound a polynomial P{z) from below on |2;| ^ 1 by the minimum of the 
cubic polynomial P mod 0{z'^) minus the sum of the absolute value of the remaining 
coefficients. Likewise, to obtain an upper bound, we take the maximum of the cubic 
polynomial P mod O(-z^) plus the sum of the absolute value of the remaining coefficients. 

Proof of Lemma \2.4\ For property (i) we first note by inspection that p'i{r) has the fol- 
lowing property: the coefficient of its first term (r term) is bigger than 1 while the sum 
of absolute values of the remaining coefficients is less than 1, implying that p'i{r) > for 
r € [0,r3]. Since obviously pi(0) > 0, we see that l/pi{r) is decreasing and positive for 
r G [0,r3]. All coefficients of P2(^3 + -2) are positive. Thus l/p2{r) is decreasing for r ^ r^. 
In particular Q{r) ^ Q(0) < 22/5. 

Property (ii) is immediate by construction. Indeed, note that pT|) defines a. Q £ Cp in 
such a way that the values of the function and the value of its first derivative match up 
at r^^r-j. The vanishing Q'{ti) = is a consequence of the fact that pi(r) has no linear 
component {a\ = 0, see Table [Tj). 

For (iii) we start with the interval [0,r2], where we will show 

(l+r)e~2>i(r) > 1/5 
Instead of showing this directly, we first notice that 

e-^^(pi(0) + ^ + r^)>l/5 
by explicitly finding the maximum via differentiation. Thus it is sufficient to show 



or equivalently 



YJr 
(l + r)pi(r)>pi(0) + — +r2 



17 

PiiC'^) := ((1 + r)pi{r) - pi(0))/r - — - r > 



For this purpose we use the partition tt = (0, r2) and Note 12.51 to re-expand pn in 
z. The coefficients of z^ , j > 3, are all very small and a direct calculation shows that 
pu > 1/50 > on [0, r2]. Therefore Q{r) < 5(1 + r)e~'^''' for r G [0,r2). 

Similarly, one can show e~^^'^pi{r) is increasing on [r2,r3] by using the partition tt = 
(r2,r3) and estimating e'^'''^(e~'^'''^pi)'(r) using Note [231 

In the interval [r3,r7] we consider the polynomial i52(r) := e^^'^{e~^^'^p2{r))". Using 
the partition tv = (r3,r7) we obtain P2{r) < —1/5 < and thus e~^^'^p2{r) is concave. 
This, combined with the fact that e~^'''='^p2(^fc) > 2/11 for k = 3,7, shows that Q{r) < 

For r ^ ry we have using Q 

(re7p3(r))' = Br-^e-^'' - 4Be2^Ei(4r) > 



10 
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Therefore re^ j-p'^ir) is increasing. Since obviously limf,_>.oo ?^e''/p3(r) = A < Hg and 
ttc^"^ /p^{ri) > ^, the estimate follows. 

For (iv), we first introduce the notation /i := [0,r3] and I2 '■= [r3,rj]. For j = 1,2 we 
let yj{r) = l/pj{r) and define 

Rjir) = -y'-{r) - 2r-^y'j{r) + yj{r) - y^{r), r £ Ij 

and consider the polynomial Mj{r) in Q[r] given by 



Mj{r) ■.= rJ-^p]{r)Rj{r) 



We introduce the partitions tti 
and define 



J_ 1 A 2 3 4 _9_\ J — I'A 11 13 83 21 5^ 

10' 5' 10' 5' 5' 5' 10^ ^^ 2 1 10' 10' 10' 50' 10' 2' 



Mj-fe(z) := Mj{ijk{z))] Pjk{z) ■■= Pj{^jk{z)); j = 1,2, \z\ ^ 1. 



(24) 



(where £jfc is the /cth component of ^j). We proceed to estimate Mjk{z) on the unit disk 
as described in Note 12.51 

This yields, for each 1 ^ A; ^ 7 and for any — 1 ^ z ^ 1, 



Mifc(z) 



Pik{zf 
On the interval [rs , r^ we have 



^ sup|Mifc(z)| ^ /ll \ /ll 

^ pi.(-i)3 < ^^ U " ^'^ V ^ '^ U " ^'^ ' 



M2l(z) 



10sup|M2i(2:)| 1 fll „ , , 



^lfe(^)P2l(2)^ 

On the interval [r4, ry] we have for each /c = 2, ..., 6 

M2fc(z) 



sup|M2.(.)| / 13 _ A / 13 _ ,,,(, 

4(-i)P2.(-i)' ^ U 2u;y-/>2^5 2u, 



4fc(2)P2fc(^;)^ 

These give the desired estimates of |-R(r)| for r € [0, ry]. 

Finally, for r £ I^ := [r7,oo) we write, with A and g{r) as in (J19p . 



(25) 



(26) 



(27) 



2/3(0 : = 



^ = ys{r; A) + {B - A^)g{r) = ^+ A^g{r) + {B - A')g{r) 
Pzyr) re 



Then the error R{r) has the form, with e := B — J^ , 

R{r) = Ro{r) + eRi{r) + e'^R2{r; e) 
where 

Mr) ■■= -y'iir- A) - -y'sir; A) + y^ir; A) - y^ir; Af 



Ri{r) := irg{r))" + g{r) - 3y3{r;Afg{r) = ^ 

R2ir; e) := -?>Ag\r) (^ + A^g{r) + ^ 



^g[r)[A—+A^g{r] 



(28) 
(29) 
(30) 



First, using Lemma 12.11 we obtain h{r-j) ^9-10 ^ (with h defined in that lemma) and 
thus 



\Mr) 



A^e 



3r 



-h{r) 



l-A^h{r) + -A^h^{r) 



-3r 



^ 3^5/1(1 



(31) 



By (fTOj) . Q > g > —e ^'"/(8r'^) and thus, noting that —r ^{rg)" + g = e ^^ /r^^ we get 

-^'^ ^ 3^2g-2rx iie-3'- 



|i?i(r)| < 



1 + 



8r2 



< 



10 r3 



Vr ^ ry 



(32) 
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Finally, 

7.3 g4^3 \^ ^ 8^3 8Ar3y 2 r-^ 

Combining ([28]), dMD and dMD with ([ISD yields 



\R{r)\ < (3A^h{5/2) + lle/10 + 3 • 10"^eV2) ^ < ^ ^ 
for all r ^ ry as stated. D 

3. The exact soliton 

3.1. Finding the exact ground state Q near the approximate one Q. We now 

need to show that 1 — Q/Q\ is small. Eq. ([3]) implies 

-§»(^r)-2^ + {l-3Q\r))6{r) = -R{r) + 3Q{r)5\r) + 6^{r) {5 := Q - Q) (34) 
r 

The boundary conditions are S'{0) = 0, 5{oo) = 0. We shall describe, again via polyno- 
mials, how to find an approximate fundamental system for (j34p . The challenge here is of 
course that we cannot hope to find an exact fundamental system for this equation, but 
require something close to it in order to set up a contraction for 6. We find it technically 
convenient to find an exact fundamental system for a homogeneous equation which is 
slightly different from the one in ([3i 



3.1.1. An approximate Green's function. 

Definition 3.1. We define ipi{r),{p2{r) as follows. Set Ji := [0, ri), J2 := [ri,r4), J3 := 
[^"4)^7)) Ji '■= ^75 00) and let 

"Pjir) ■■= 9i+2(3r)xoi + qj+2{r - l/2)Xi4 + qj+2{r - 2)^47 + aje^'^'Xioo (35) 

for j = 1,2, where ai = 1,02 = 1/2,0"! = —^,(^2 = 1, Qs^Qi are of the form qj{r) = 
X]fc=o He ' where the 6^ are given in Table [71 The factor ^ in front of e'' in the definition 
of ip2 is chosen so as to normalize a Wronskian to 1. Finally, set Qj = fj on J4 and 

9j(.r) := ^j{r) + gj{r',+) - if,{r',-) + {g'^{r',+) - 'f'j{r',-)){r - r\) V r G J,_i (36) 

for all i = 2,3,4, where J^_i = [r'^_i,r'g) (we use this notation only for ([36]) ). 

We remark that the jumps appearing in (|36p are very small; more precisely, 

max{\ipj{re+) - ^j{re-)\,\ip'j{re+) - ip'j{re-)\) < 3 • 10"^ V £ = 2, 3, 4, j = 1, 2 

The functions gi,g2 from the previous definition satisfy an ODE which is a perturbation 
of our main Sturm-Liouville equation. 

Lemma 3.2. The functions 51, (72 are in C"'^([0,oo)) n Cp, and solve the ODE 

-y"{r) + U{r)y'{r) + V{r)y{r) = (37) 

where U and V are piecewise rational functions which obey the estimates 

||f/||oo < C/'^ := T^, \\V-l + 3Q^\\oo<V':=^ (38) 

loo 3d 

as well as U = 0, V = 1 on J4. 
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Proof. That gj are in C-^([0,oo)) n Cp is clear by construction. One has for r £ Ji 
1 ^ ^ ^ 4, 

jj^^ ^ 9'i{r)g2{r) - g2{r)gi{r) ^ _ g'i{r)g'2{r) - g2{r)g[{r) 



g'l {r)g2 (r) - g'2 ir)gi (r ) ' g[ {r)g2 (r) - g'^ (r)gi (r) 

which are manifestly rational functions. Moreover, [/ = 0, y = 1 on J4. To obtain the 
bounds of dMD we rely on Note [231 Let vr = (O, ^, ^, ^, |, f , f ' 1' i> 2, f)- 

To estimate U , one writes within each subinterval U{r) = M(z)/P{z) with polynomials 
M,P G Q[e'^^e-''7][z] and obtain 

\U[r)\ ^ ^(^L 

' ^ ^' |P(0)|-P(1) 

where M[z) is obtained by placing absolute values on all coefficients of M, whereas P is 
the result of the same procedure applied to P{z) — P{Q). For V , one proceeds similarly, by 
first writing Pj{r)'^{y — 1 + 3(5^) as a rational function on (0, r-j) with j = 1, 2 depending 
on whether r < rs, or r > r-^, respectively. Substituting the same affine change of variable 
as for U into the numerator and denominator of this rational function now establishes the 
desired bound. D 

We shall need to modify the system gj to accommodate the -^ term in the three- 
dimensional radial Laplacian. In the following lemma, note that g^ is regular at r = 0, 
but grows exponentially, whereas g^o decays exponentially, but is singular at r = 0. 

Lemma 3.3. The functions 

go{r) := r~^ ig2{r) - ^7^9i(^) j ' 9oo{r) := r~^gi{r) 
form a fundamental system for the equation 

- y"{r) + [u{r) - -^ y'{r) + [v{r) + ^\ y{r) = 0, 

and their Wronskian satisfies the estimate 

W{r) = r^^ if r '^ r-j 

< r^\W{r)\ ^1 V < r < ry 



(39) 



(40) 
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One has the following pointwise bounds: 

\goo{r)\ ^ (^ + ^)Xo4(0 + ^Xiooir) and \gUr)\ ^ ^Xo^ir) + ^t 
as well as 

\go{r)\ < ^Xjg i^(r) + l^X^i^^^ir) and \gQ{r)\ < 18xo4(^) + ffX4oo(^) 

Proof. The first claim follows immediately from ()37p . Let gj{r) := r^^gj{r) for j = 1,2. 
Then the Wronskian W{r) = go{r)g'ao{i') — g'Q{i')goo{r) of ^oo and 50 satisfies 

W{r) := ~gi{ryg'2{r) - ~g[{r)g2{r) = r-^ e^n""^'^"' (41) 

This follows from the fact that W is continuous by Lemma 13.21 and 

W'{r) = {U{r)-2/r)W{r), W{r)=r~^ \i r ^ r7 

by dSSD. From ([MD one now obtains (|iO|) . 

For the estimates of g^ and goo we use re-expansions and the following simple observa- 
tion. 
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Note 3.4. Estimates of functions / for which /' is a quadratic polynomial multiplied by 
a monomial of any degree or by an exponential are elementary. We can use this in our 
estimates as follows: if / = YlT=i fj ^^^ J C M then, clearly, 

m m 

inf/^^inf/,; sup/^^sup/,; (42) 

3=1 ^ i=i ^ 

If / is a polynomial we write / as a sum of subpolynomials /j, the first one containing the 
monomials of degree ^ 3 and the others consist of the monomials of degrees G [3Z+1, 3Z+3] 
for all / with 3/ + 1 ^deg/. If /(r) = e P(r) , the same applies to a decomposition e Pj 
in which Pi has degree two and the Pj , j ^ 2 consist of the monomials of degrees 2j + 1 
and 2j + 2. 

To establish \gi{r)\ = \rgoo{r)\ < | + 1^ on (0, 1) let mi(r) = l + %- The result follows 
using Note 13.41 for the polynomials gi it mi taking ^ = ^ + I5 and re-expanding (as a 
polynomial in s) s G [—1, 1] and then r = | + s, s G [— ^, §]• 

To establish |(7Jxd('^)I < ^ o'^ (0, ?"4), one proves the equivalent 

\rg[{r) - gi{r)\ < - < r < r4 (43) 

The result follows once more from Note [331 re-expanding rg[{r) — gi{r) via ^ = Jj + |g 
for \s\ <l, r = 1 + s, s £ [--g, ^], and r = ^ + s, s G [0, ^]. 

Next, we turn to the estimate 1500(^)1 < - — on r ^ r4. On r ^ rj we have an exact 
equality to - — so it suffices to deal with r^ ^ r < rj. We note that the inequality is sharp, 
and we need a different method: we let ip = e^gi and look at ip" . Here we use Note [33] 
and expansions at rg + s, |s| ^ j^, and r = 2 + s, s G [— ^, ri] and r = rs + s, s G [— ^, n], 
to see that ip" < 0. Since (p'^r^) > and ^p'irj) = 0, we see that tp' ^ 0. Since ip{rj) = 1, 
the property follows. 

Finally, one has the bound |5'^(?')| ^ (1 + r)e~''/r^ for r > r^. In view of the exact 
expression one has for r ^ ry it suffices to deal with r^ ^ r < r^- Thus, we need to verify 
that 

e^\rgi(r) — gi(r)\ ^ 1 + r, r^ ^ r < rj 

We let (/Ji(r) := e^'{rg[{r) — gi{r)). Using the partition tt = (r4, 2, 5, jQjTy), explicitly 
minimizing the leading cubic polynomial and taking the £^ norm of the rest, we see that 
(Pi{r) > 0. Thus, by monotonicity, (Pi{r) ^ v'i{'''7) < — 1- This implies that (pi{r) < 
on [r4,r7]. Moreover, {pi{r) + 1 + r ^ on [r4,r7] since ipi{r) + r is decreasing and 
ipirj) + ry + 1 ^ 0. 

For go one proceeds in a similar fashion. We begin with the bound 150(^^)1 < 13/2 on 
(0, |). On [0,ri] we apply Note 13.4] for / with r = -^ + ^s while on [ri,r4] we look at 
the polynomials rgo{r) ± ^ with r = ^ + s and ■? G [— ^, ^]. 

Next, one verifies that 1(70(^)1 < |i+f on r ^ 2- 0^^ the interval r ^ r-j one checks the 
explicit expression gQ{r) = %- + —r where k = — f^oj G (0,4) and thus 



2r "•" re^ wiicic n. -^^j^ 



. > 6 K 

9o{r) = TT + 



2r re'' 1 + r 



1 + r /I k 

2^i2F 



7/1 k \ e"" 7 



5V2 100M+r51 + r 



On the interval [21^7] 0116 can check via the partition tt = (2)l5|)i) that g^ > 0. 
Furthermore, we apply the same partition to the expression 

5(52(0 +%i(r))(l+r)-7re'' 

which is of the form admitted by Note 13. 4i One then sees that it is negative. 
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For the bound |(7o('")l < 1^ o^ (Oj'^i)) we multiply through by r^. We then use the 
substitution r = ^ + s. On [ri, r4] one uses r = ■^ + s. 

Finally, we verify |5o(^)l < ff for r ^ r4. On r ^ rj one checks that 

-r / / \ 1 -2r -'^ Z-' -2r 

re qn(r) = ae e 

with a,/3 > 0. The right-hand side is clearly increasing in r and < 2- Since goirj) > 0, 
the claim holds for r ^ ry. On r4 ^ r < ry we re-expand rgQ{r)e~^ via r = | -|- s with 

Now we come to the main result of this section, which is the estimate of the relative 
error between Q and Q. 

Proposition 3.5. Let Q be the exact ground state of dSD and Q be the approximate one 
given in Definition \2.3[ Then one has the error bound 

\Q{r) - Q{r)\ ^ eo^-^ V r ^ 0; eo := 7 • 10^^ (44) 

1 + r 



-5"{r) + (U{r) - ^y'{r) + (y(r) + ^)<^(r) = /ii(<5,r) (45) 



Proof. Rewrite ()34p in the form 

'-)6'ir) + (vir)+'^^^ 
r J \ r 

where 

/i^(5, r) := -R{r) + U(r)b'(r) + \v{r) - 1 + -iQ^ir) + ^^ b{r) + 3Q(r)(5^(r) + b\r) 

V r 



We seek a solution to (j45p which obeys the boundary conditions 

<5(0-F) G M, (5(oo) = 
In fact, this solution is unique and is of the form b = H{6) where 

H{6){r)=g^{r)J^ ^^^^ ds + g,{r) J^ ^^^^ ds 

=: H,{6){r) + H2{6){r) 
in terms of the fundamental system from Lemma 13. 3[ We also have 

gois)hi{6,s) , f°° goo{s)hi{6,s) 



(46) 



[H{6){r)\ =: H = g^{r) J — ds + go{r) 



■ ds 



W{s) 
=:H[{6){r)+H',{6){r) (47) 

Lemma 3.6. H is a contraction in the ball 

X:={fGC\iO,oo))\\\f\\x<eo} (48) 

with norm 

||/|U:=sup(r + l)e^(|/(r)| + i|/'(r)|) 

Thus there is a unique fixed point Sq £ X . 

Corollary 3.7. Therefore, y '■= Q + Sq > solves ([3]) on (0, oo]j, remains bounded as 
r —7- 0-|-, and decays as r —t- oo. By Coffman's theorem |10] . this uniquely characterizes Q 
whence Q — Q = 5q. 



Since y €z C^ hy construction, it is a weak solution of ((3} and by standard Sturm-Liouville theory 
therefore also a smooth one. 
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3.2. Proof of Lemma 13.61 For any r G [OjTy], denoting a;(r) := e "^ /{I + r), any 6 G X 
satisfies 

\hi{6,r)\ ^ \R{r)\ + ||C/||oo|5'(r)| + {\\V - 1 + SQ^^o + r-i||C/||oo)|5(r)| 

+ mQimioo + \\s\\imr)\ 



^ \Rir)\ + 



eo 



1 



1 



33 V 36 165 ' ° ° 



uj{r) ^ |i?(r)| + 



3eo £o 
50 165r 



(49) 
and [231 with 



For r > ry we state a different bound, see the case r > ry in Lemmas 
„ ._ _9_ 
^3 ■ 200' 

|/ii((^,r)| i^^ + 3Q\r)\5{r)\+3Q{r)6\r) + \6Hr)\ ^(^3 + 24^0)^ V r > ry (50) 

We now show that H takes X to itself, which is based on the estimates of the fohowing 
lemma. This is the most technical part of the contraction argument for 6, and it is proved 
by inserting ([32D and <^U^ into ([MD and using the estimates from Lemma [331 We suppress 
the argument r for the most part in the formulation of the following lemma. 

Lemma 3.8. Let 5 G X be fixed and consider (|46p . For all r > one has the estimates 
\H^i6)\ < (ii . 10- V^ + §) Xp 1) + (|j - ir) (8 • 10-6 + f) x^i^,^ 



+ Hi • 10-^ (3 - r) + iig) X47 + ^ (10-^ (i, - 36e--) + 



Tea 
50 



48£o 
25 



-2rN 



X700 (51) 



|i7l(<5)| < [10-V (2 - f ) + §] X,, K + (i - 3r) (8 • 10"^ + §) x,i 



[04) 



,1) 



+ HI • 10-6 (3 - r) + ^eo) X47 



(r-+l)e- 



:iO-Mi,-36e--) + (^-ie--).o)x7oo (52) 



Hm < [10-^ (i - y) + §] X. K + f§ [10-6 (i + 2 (1 - r)) + ^.0] X,l , 



[ii) 



+ 5(;^ (^ • 10-^ + 10-^(6 - ^0^0) X4r + ^^S^gl^Xroo (53) 



F^(5)| < ft [10-^ (i - y) + f|] X,, U + 18 (10-6 (ii + 2 (1 - r)) + ^.0) X,i 



loi) 



1) 



+ 



r+l UO 



. 10-6 + 10-2(6 - ^r)eo) Xa7 + 



3e-3^(3+1600£o) 
1600r^ 



X700 (54) 



Proof. To avoid working with Ei(x), we write co (s) ^ e^/(a + 1) inside every definite 
integral from a to b, and e~^^/s^ ^ 4e-^*/(25s) for s ^ rj. Also, estimating sums or 
products is of course elementary: for instance, e^ + ax^ + bx + c < d is equivalent to 



-^(d 



ax' 



+ 5x + c)) — 1>0- checked by examining the derivative. Using this, the 



result follows by straightforward calculations of the integrals using Lemma [37 
([50|) . See Section [6l for details. 



9|), and 

D 

Now we can show that H (5) takes the eo-ball of X to itself 

Corollary 3.9. Let {Ho{6)){r) = \H (6) (r) | + l\H' (5) (r) |. We have 

(r + l)e'"(Fo(5))(r) s^ eo Vr G M+ (55) 

Proof. This is a straightforward consequence of Lemma 13. 8| the details are given in ^6.31 

D 
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3.3. Contractivity of the map. Let (^o(^) = "^il^) ~ ^2{f)- Clearly, 5i^2 only occur in 
the integrands in Ho{6i{r)) — HQ{62{r)) (through hi{6i,r) — hi{62,r)). Note now that 



\hi{5i,r) - hi{62,r)\ ^ {\\V - 1 + 3Q''\\^ + r-'\\U\\oo)Mr)\ 

+ (6||Q||||5o||oo + 3Po||L)l'^o(r)| + \\U\\ooK{r)\ ^(^ + 



1 



50 165r 



oj{r)\\6o\\ (56) 



If we replace eo by ||5o|| ^ind set R = 0, the last term in (|49p is identical to the last term 
in (j56p . Hence, with the replacements above, all contractivity calculations shadow those 
for the bounds on hi. Thus, with virtually the same proof as that of Corollarv 13.91 see 
Section 16.31 one derives the estimate 



(r + l)e^{\H{6i){r) - H{62){r)\ + \H'{5i){r) - H'{62){r)\/5)/\M 



^ 4e^(r + l) / 13 llr\ , ^^ .^ 

^ ^2^^[o,i) + ilOO - ^J (" + '^' ^[^D + 



lO^^gr /37 ^ ^2r'^ 



+ 



21 7 

+ 



49 _2,. 



r 
431 



h9r 

20 



67 \ 



X47 



23 3 

10 "^ 5^"^ 50r2 25r7 



125 250r^ ' 250r ' ^ ' in ' k„^ ' c;n»,2 ot;». / ^7oo / V / 

where the final bounds follow by differentiating in r for the first three and in 1/r for the 
last. For the last one we note that — f§ + ^ + ^^ — ^ is negative and decreasing. D 

3.4. Further estimates. In the study of L± we need a sharper estimate oi Q — Q. 

Lemma 3.10. For r ^ -^ = r^, we have 

6{r) = Q{r) - Q{r) = hie^' /r + b2{r) 

where \bi\ < 5 • 10^^ and \b2{r)\ < ;|j^;:t-- 

Proof. It follows directly from (|46p that 



6ir)=9Ur) I ^M^^^ + ^oolO f ^A^M^ ds 



+ 9o{r) 



W{s) 
goo{s)hi{6, s) 
W{s) 



W{s) 



(58) 



ds 



It then follows from (j50p that 
'■'=^go{s)hi{S,s 

'0 

and 



W{s) 



■ ds 



9oir) 



10"^ 


9ooir) / 


goo{s)hi{6, s) ^^ 


W{ 


s) 



go{s)hi{6, .s) 
Wis) ^^ 

1 e-3^ 



< 



1 e- 



-3r 



25 H 



50 r3 



(59) 

(60) 
D 



which concludes the proof. 

4. The operator L+ 

In this section we prove the first of our main results, namely the gap property of L+. 

Theorem 4.1. The operator L+ has no (L'^{M.~^)) eigenvalue or resonance in [0, 1]. 

Standard ODE analysis shows that there are two solutions ui{r;X) and U2{r;X) of ^ 
with the properties ni(0; A) = 1 and U2{r; A) = r~^e~^^^~'^(l + o(l)), r — )• oo. These are, 
up to constants, the only ones acceptable at zero and infinity respectively, see §6.11 

Let W = uiu'2 — U2u'i be the Wronskian of the two special solutions ui^2- As mentioned 
in Section 11.31 the existence of an eigenvalue or resonance of -L+ is equivalent to W = 
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for some A. Note that W is not constant due to the first order derivative in L±. However, 
W'{r) = —-W{r) whence VF(ro) = at one point imphes that W{r) = everywhere. 
Theorem 14. II is a corollary of the following result. 

Proposition 4.2. We have the following estimate 

inf \W(rT,X)\ 5? 43 • 10"^ (61) 

Ae[o,i] 

4.1. Proofs. As mentioned in Section [L3l we construct the quasi-solution wf, a piecewise 
polynomial of two variables r and A, on the interval [0, ry], and Wg whose expression 
involves exponential integrals on the interval [r7,oo). To show that wf is close to ui 
and W2 is close to U2, we use the same contractive mapping strategy we used for Q — Q 
outlined in Section 11.31 the only difference is that the equations and solutions depend 
on the parameter A. The method for obtaining the estimates is explained in Note 14.31 
below. Finally, to estimate the Wronskian, we first approximate tf;^(r7) and drW2{r7) by 
polynomials of degree 7 in A, and the calculation then reduces to estimating a polynomial 
using Note 12.51 

4.1.1. A polynomial quasi-solution for r ^ rj. To build the quasi-solution, we first define 

wi{r)= Yl XjCki;jX''z^ (62) 

(j,fc,/)6S 

where S = {j,k,l \ 1 ^ j ^ 3, ^ k ^ Mj, ^ / ^ -^j}) (^kij are given in the appendix, 
Xj, j = 1,2,3, are the characteristic functions of [0,r2), [r2,r5), and [r5,r7], respectively, 
Mj ^ 15 and z depends on r as specified in the top rows of the tables in the appendix. 
For the most part, we suppress the A-dependence in our notations. To ensure that wf is 
C^ we next let 

wf{r) = wi{r)x57 + [wi{r) + Wiir^) - wi{r5-) + (u)i(^5) - w[{r5-)){r - r^)] X25 



+ 



wi{r) + 'wtir2)-wiir2-) + ^{drwt{r2)-w[{r2-)){r^-4) Xo2 (63) 



Note 4.3. In this section, we divide [a, b] x [0, 1] "vertically", using a partition tv of [a, b] 
in r and the trivial partition (0, 1) in A. Clearly, such partitions are determined by tt. 

(i) Let P{y,z) be a polynomial with {y,z) G [—1,1]^ and Pi defined to be the sub- 
polynomial consisting of all terms of P of the form CkfiZ , ^ A; ^ 3 and cojy^ , 1 ^ j ^ 3. 
We bound P above (below) by the maximum (minimum, resp.) of Pi plus (minus, resp.) 
the sum of the absolute value of the coefficients of P — Pi. Since Pi is the sum of a 
cubic polynomial in z and a cubic polynomial in y, the extrema calculations reduce to 
solving one variable quadratic equations. We proceed in this way for simplicity, given 
that the mixed terms typically do not contribute too significantly to the final estimates 
of our polynomials, (ii) If /(r) = e'^^Q{r) (equivalently, /(r) = e"'^ /Q{r)) where Q is a 
quadratic polynomial, then /' = iff a quadratic is zero; thus e"^ it Q(r) > c where c is 
any constant can be verified by elementary means. 

Lemma 4.4. We have \wf{r)\ ^ M{r) := — ]^q~ X02 + ^^27 

Proof We introduce the partition of [0, r-r] x [0, 1] induced by tt = (0, 3/10, 17/25, 3/2, 5/2). 
The lemma now follows by applying the method in Note 14.31 to wf{ik{z), 2(1 + y)) i 
M{ik{z)) for 1 s$ /c ^ 4. D 
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Now we invoke the estimates for |Q — Ql (Proposition I3.5|) and Q (Lemma l2.4|) to obtain 
|3Q2(r) - 3g2(r)| ^ 3(2Q + \Q - Q\)\Q - Q\{r) 



^ 1 



10000 



21e-=^^X02 + ^(l + O-'e-i3-/5x27 



+ 1 . iQ-\r^r 



< 1 



5000 



2 

lle-3^Xo2(0 + 12(1 + r)-ie-i3-/5^27(r)] (64) 
We shall use this estimate repeatedly, for example in the following bound on the remainder 

R^ ■= -w+" - 2w+'/r + (1 - A - 3Q2)y;+ (65) 

Lemma 4.5. We have 

\Ri{r)\ ^ 10-^ {§(11 - 35r + Mt^)x,2 + ^X25 + ^Xsr} (66) 

Proof. We first consider Rq = —wf" — 2wf /r + (1 — A — 3Q^)wf and the piecewise 
polynomial 

Roir) = [Ro{r)/Q\r)]xo2 + [rRoir)/Q\r)]x27 (67) 

On [0,r2] we use the partition tv = (O, j^q, 3^, 20) 25' Too' 50' ^2)- Since Q is positive 
and decreasing, in each subinterval we only need to show that \RQ{ik{z))\Q^ {ik{—l)) is 
bounded by the right-hand side of (f66|) . It is easy to see that this can be reduced to 
applying the method in Note SSI for ^o(4(-z))'3^(4(-l)) ± 10-^(11 - 43r + 54r2). In 
this fashion, we obtain the estimate 

\Romz))\ ^ \Romz))\Q\ik{-l)) < 10-4(11 -43r + 54r2) 

for all 1 ^ A; ^ 7. 

Next we introduce the partition of [r2,r5] provided by tt = (^"2, fg, '"Sj fg, fgj'^s) to 
estimate |^o(4(2))|Q^(4(-l)) ± 2 • 10"^4(z). This yields the estimate 

4(^)|iio(4(^))| ^ |i2o(4(^))|Q'(4(-i)) < 2 • w-'h{z) V 1 ^ A: s; 5 

For the interval [r5,r7] we use the partition induced by tt = (^5; 3^) x§) ff; ffj ffj fg;'"?) 
and concludes that 

ek{z)\Romz))\ ^ |iJo(4(^))|Q'(4(-l)) < 15- lO-64(^) 

We have |-Ri — Ro\ = 3\Q^ — Q'^\ \wf\. Combining these results with (I64p and Lemma 14.41 
to estimate \3Q^ — 3Q'^\\wf\, we obtain 

\R^{r)\ ^ 10-4 { [11 - 43r + 5Ar' + Hiilzlle-^'^] ^02 + (l + ^^^rS^) ^ 

+ (ll + fS^)x57} (68) 

The right-hand side of ()66p minus the right-hand side of (|68p is positive; this follows using 
Note lTSl (ii), after the substitution of 5e"'^''xo2 by the bound (5 — llr + 7r^)xo2- ^ 

4.1.2. A quasi-fundamental system of solutions on [0,r7]. To show there is an actual 
solution ui of (JH on [0, ry] with 6 := ui — wf small, we construct two functions gf2j 
approximating two linearly independent solutions of ([!]); define first 

9tir)--= E dt,,X,X'z'; ~9tir)--= E 4i;X,X'z' (69) 

(i,fc,0e5 {j,k,i)es 

where S = {{j,k,l) | 1 s: j s: 3, s^ A; ^ Mj, s: / ^ 15} and now Xj, J = 1,2,3, are 
the characteristic functions of [0,ri), [ri,r4), and [r4,ri], respectively. The coefficients. 
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the intervals corresponding to j = 1,2,3, and the expressions z = z{r) are given in the 
appendix. 

To ensure C^ behavior we use the same method as in Definition 13.11 to set g^ = g^ on 
Ji and 

5+(r) := 5+(r) + ?+(r,-) - ~gf{r,+) + (ff+'(r,-) - ~gj\r,+)){r - n) ^ r e J, (70) 

for ah I = 2,3,4 ancO 9t{f) = dtif)/^- In contrast to Definition 13.11 where the C^- 
matching is done from right to left, we find it necessary to carry out the matching by 
going from left to right. 

We construct a second order equation satisfied by gf, 5^ in the form 

- g" + {A{r) - 2/r)g' + {B{r) + A{r)/r)g = (71) 

where 

^(^) ^ 9tir)U"ir)-Uir)9t"ir) ^^^ ^^^^ ^ gl'{r)U"{r) -U'ir)9t"ir) 
gtir)gf'{r) - gf{r)g^'{r) g^{r)g^'{r) -gf{r)g^'{r) 

Lemma 4.6. We have \A{r)\ ^ 1^(5X06 + 4X6?) and \B{r) - 1 + A + 3Q^{r)\ < ^. 

Proof. We use a partition tv = (O, ■^, ^, ri, ^, ^, ^, |§, "^4, |, t§, ^-,^6, n) and estimate 
above/below, by re-expansion, the absolute value of the numerator (denominator, resp.). 

□ 

We rewrite the equation for 6 (see beginning of §4.1.2p in the form 

- S" + {A{r) - 2/r)6' + {B{r) + A{r)/r)6 = Ri{r) + A{r)6'{r) 

+ [B{r) - 1 + A + 3Q{rf + A{r)/r + 3{Q{r)^ - Q{rf)]5 =: /i2 
Note. In the following we write ||/|| for suprg,. 1 |/|. 
Lemma 4.7. We have 

1/12(5, r)| ^ 10"^[33(ll-35r + 34r2)/10 + 5||5'|| + (82 + 5/r)||(5||] Xo2 
+ 10-4 [(11 - 6r)/5 + 70PII + 5p'||] X25 + lO^^ [(39 - 9r)/100 + 64||5|| + 5||5'||] ^56 

+ 10"4 [(39 - 9r)/100 + 78||5|| + 40||5'||] Xer (72) 

Proof. This is obtained by combining (j66p . Lemma [4.6l and (j64p and using the monotonic- 
ity of the coefficients containing exponentials. D 



11 
10 



Lemma 4.8. We have \gf\ < i^Xgi + ^Xi7 and bjK ^ + 

Proof. We use the partition induced by vr = (0, ri, ^,r4,2,r'j) and apply the method in 
NotelUto Xoi(Tn ± at), Xi7(i ± ^at) and il + ^ ± rg+. D 



In addition, one has the following. 



Lemma 4.9. There are the bounds \gt'ir)\ ^ 3xo2 + X25 + ^^5?; bj'(^)l ^ Ti^Xo2 + 



f2-X25 + T0^57- 

Proof Similar, using tv = (0, n, r2, r4, rs, j^, fg, ry). D 



We note that g^ is approximately ui]^ redefined on the same interval as gj for convenience. 
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4.1.3. The actual smooth solution on [0,r7]. Let 

u (K\ +( \ r 9t{s)h2{5,s) + r g^{s)h2{6,s) 

Hoi5) = gtir) I ^(,) ds - atir) j^ ^(,) ds (73) 

Clearly, we have 5 = Hq{6) and 5' = Hq{5) where 

HoiS) = 92 (r) J^ ^(^^ ds - 9, (r) j^ ^(^^ ds (74) 

Lemma 4.10. There is the bound 

\Hoi6){r)\ + \H',{5){r)\/5 ^ 1/1200 + ||5||/10 + ||y||/80 (75) 

Thus H is a contraction in the ball 

X := {/ G C\{0,rr)) \ \\f\\x ^ 1/1080} (76) 

where \\f\\x ■= sup.^g[o^.^^) (|/(r)| + \\f{r)\) (cf. also footnote \3l\ on p. [7^ . 

Proof. We crudely estimate the quantities Hq and Hq by placing absolute values on all 
terms, and by using the bounds already calculated for 9^2 etc; \6\ and \6'\ are estimated 
by their supremum norms, written as || • ||. Since W{r) := gi~(?')(7^'(r) — gi'{r)g2{''') = 



r 



have 



exp ( Jq A{s)ds j ^ r ^ exp ( — f^ \A{s)\ds j (we note that W{r)r'^ — )• 1 as r — ;• 0) we 



l/|VF(r)| < -TTT^r"' < 



625 50 

On the first interval, [0, ri) a direct calculation shows that HQ{6){r) := \HQ{6)\+\HQ{5y {r)\/5 
is majorized by 

max {\rP^{r)\ + \r P2ir)\\\6' \\ + \P3{r)\\\5\\) ^ 2 ■ 10^^(200 + 1300\\5\\ + 57||5'||) (77) 

where Pj are polynomials of degree j, easily maximized since they are increasing on this 
interval (all have positive coefficients except rP4; (rP4)' has positive coefficients after we 
replace r^ by rri and r^ by rrf); thus rP^ is increasing. 

Calculating HQ{6){r) for r € [ri,r2), we obtain a rational function; we first replace 
1/r, 1/r^ by 1/ri, 1/rf respectively (their coefficients are positive) and we get an expres- 
sion of the form P5(r) + P3(r)||(5|| +P3(r)||(5'||, with the same convention as above for the 
polynomials (different from the Ps on the previous interval). Once more, all polynomials 
except P5 have positive coefficients. In P5 we first replace r^ (whose coefficient is positive) 
by r^r2', the derivative of the new polynomial has an explicit positive minimum. Thus 
the maximum of P5 is reached at r = r2. Thus, by taking r = r2, we get, for r G [ri,r2), 

^„(,)M<^^ + M + M (78) 

^ ^^ ^ 25000 125 2500 ^ ^ 

On the interval [r2, rs) we proceed in the same way, replacing 1/r, 1/r^ by l/r2, l/r| resp. 
This results in an expression of the form P4 + P3||(5|| + P3||(5'|| with the same properties and 
conventions as above. Now the derivative of P4 can be minimized explicitly: it is positive 
and thus P4 is maximal at the right-hand endpoint. We get the following majorization of 

l^o(<^)(r)| 

^o(<5)(r) ^^ + — + — VrG[r9,r5) (79) 

^ ^^ ^ 1200 34 515 ^ ^ ^ ^ 

On the interval [r^jr^) we replace 1/r by l/rs and obtain an expression very similar to 

the one on [r2,r^). It is dealt with in the same way, whence 

Hoi5)ir) < 1/1250 + ||<5||/11 + ||5'||/164 V r G [rg, rg) (80) 
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Finally, on [rQ,rj), after replacing 1/r by 1/rg in the positive terms and by l/ry in the 
negative ones, we reduce Ho{5){r) to the form P4 + PsH^H + P3||(5'||. P3 and P3 are 
manifestly increasing and minr^g^^^) P4 > 0. Thus the maximum is reached at r^ and we 
get 

^o(5)(r) ^1/1200 + ||(5||/10 + ||5'||/80 VrG[r6,r7) (81) 



Of all estimates, the worst bounds are in (ISip : contractivity as well as preservation of the 
ball thus follow from dH]). D 



Corollary 4.11. The function wf differs from an actual solution uq of L^u = Xu by at 
most 1/1080 in \\ • \\x- Furthermore, 

\6'{rj)\ s: 1/1800 

Proof. The first part is just ()76p . The second part comes from direct substitution in \Hq\ 
followed by Lemma 14.101 

\Ho{Sy{r7)\ ^ 1/2080 + ||(^|U/17 ^ 1/2080 + 1/17 • 1/1080 

as claimed. D 

Lemma 4.12. We have \wf{rY)\ < 19/47, |wj*"'(r7)| < 4/19. 

Proof. Viewed as quintic polynomials in A, wf{r7) and wf{rj) are simply estimated as 
in Note 14.31 (reduced here to one variable y where X = ^{1 + y)). D 

4.1.4. The quasi- solution bounded on [r7,oo). Let a = \/l — A and ai = 1 + cr. We let 
ao = 1413/64. In this region, we look for Uoc in the form W2 + Sq where 



wt{r) = ^-[l + aofiir)y, a/i(r) := -Ei(-2r) + e2-ViEi(-2air) (82) 

is close to an exponentially decaying solution of (j3|)|_| Then, u;^ satisfies 

- w+" - 2w+'/r + h - A - ao^;^] w+ = -ag ^3 /i(r) =: R2 (83) 

Lemma 4.13. We have (i) |w^(r)| < e^"""' /r and (ii) \R2{r)\ < ^^ 
Proof. By straightforward algebra, using the bounds in ([9]) one obtains 

— 2r 

f[{r) = 2aie2-^Ei(-2air) + - > (84) 

r 

Since limr_^oo /i('") = 0, we clearly have fi{ri) ^ fi{r) < 0. Now 

fi{r7;c7) = r'^lt±^du ^ f,{rr,0) = 6Ei(-5) + e"5 > --^ 
and the estimates follow. D 



To obtain this approximation, we rewrite Q as (— -rr ~ 7^ + 1 ^ '^)'^ = 3Q^u, replace Q by the 
leading term of Q, and iterate the associated integral equation. 
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4.1.5. The equation for 5q. The difference u^ 



Wn 



6q satisfies the equation 



-2r 



- 6'^ - 26yr = R2 + {X-l + ao-^)<5o + ( 3Q 
Note 4.14. One has (i) 



-2r 



ao- 



(u;++(5o)=:/i3(r) (85) 



\Q\r) - Q\t)\ < (2Q(r) + \Q{r) - g(r)|)|Q(r) - Q{r)\ ^4-10- 



-2r 



(n) 



-2r 



\3Q\r)-ao- 



^\3Q\r)-3Q\r)\ + 



ao 



—2r 

V-3Q'(0 



< 



1 e 



-2r 



20 r^ 



Indeed, (i) fohows from ^ andi^. (ii) uses ([22]) and (03]) and (i). 
Now using (j85p . Lemma 14.131 &iid Note HTT^ we have 



1^3(01 ^ l^2| 



-2r 



3Q^ - aQ- 



\w+{r)\ + \X-l + 3Q^\\6o{r)\< 



13e- 



-2r 



100r3 



+ Mr)\ 



Looking for exponentially decreasing solutions, we write (j85p in the integral form 



/°° lit f°° 
^y S^h3{s)ds 



(86) 



4.1.6. T/ie actual solution on [r7,oo). 

Lemma 4.15. Hi is contractive in the ball {/ 
{/I 11/11 =sup,>,,rVn/(r)|<oo}. 

Proof. Noting that (s/t) > 1 and t > r we write 



^ 13/300} in the Banach space 



1 



\Hi{5o)\{r) ^ -, 



dt 



'2r 



s-^h3{s)ds ^ -^^ —- + 



13 



r3 J.^ J^ r^ V400 

whence the claim. 

Hence ||-H'i(5o)|| ^ 13/400 + ||(5o||/4 and the claim follows. 
Corollary 4.16. We have 

1 3 2^ 

\u2{r,) - wt{rj)\ = \52{r-,)\ ^ ^^,e~' < 2 • 10"^ 

14(^7) -^+(^7)1 = 1^2(^7)1 < 4- 10-^ 
Proof. Only 6' needs to be estimated; this is immediate: 

K{r)\ = |(^i5o)'(r)| ^-,1 s%{s) rf^ ^ ^ (^ + ^J < 4 • 10- 
and we are done. 
4.1.7. The Wronskian of the well-behaved quasi- solutions. Eq. ([82]) implies 

26^'"^'" 



(87) 
D 



(88) 



D 



wjirj) 

2e-'""' 
25cr 



5cr 



(o- - aoEi(-5) + aoe^'^<TiEi(-5(Ti)); w'^'{rj 



(o-(2aoe"^ - 2 - Sa) + Ei(-5)ao(2 + Sex) + OQe'"' {f>cr^ + 3a- 2)Ei(-5(Ti)) (89) 
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Lemma 4.17. (i) Let z = 2a — 1. The functions wf and wf' satisfy the estimates 
79 737 z 147 z^ 49 z^ 103 z^ 131 z^ 13 z^ 6 



691 2580 2 ' 412 22 15523 ' 55524 141925 ' 340 2^ ~ 445 2^ " ^2 (^7) 

52 26 z 10 z^ 80 z^ 41 z^ 41 z^ 66 z^ 13 2^ +, 
'509"''T852~3T322~ 857 23 ^ 307 2^ ~ 392 25 "*" 1109 2^ ~ 480 2^ ^ "^^ {^7) 



<3-10-^ 

<6-10"^ 

(90) 



(ii) We have \w2{r'/)\ < 21/50 and |w'2^'(r7)| < 19/100, which will be used in estimating 
the Wronskian of the two possible eigenf unctions. 



Proof. The polynomials in ()90p are simply truncates of the Taylor series of the func- 
tions involved. These, and the estimates, are obtained as follows. Denoting w'^{a) = 
aw2{rT,a) we have ^^(ry;^) = Jq w^'{as) ds. We then approximate wf{as) using a 
Taylor polynomial around o" = ^ with rigorous bounds for the remainder, using Cauchy's 
formula. This can be obtained by expanding the exponential functions and expanding 
the integrands and then integrating the series term by term. For example, 'L^-u '^^^' 
fers from the sum of the first 12 terms of its Taylor series expansion in u at u = | 
by no more than 2 • 10"'^ (this follows from by Cauchy's integral formula). Hence 
Ei(— 5 — 5<Ts) = Ei(— 10) + L^^ ^-^-u ^^ differs from the integral of the Taylor poly- 
nomial by no more than 10~^. After obtaining a polynomial approximation of W2{r7) 
in this way, we re-expand it in z = 2o" — 1, and it so turns out that the coefficients of 
z^ for k > 7 are manifestly small. Discarding them and using rational approximations 
of the remaining coefficients we obtain the polynomial in Lemma 14.171 The result for 
u'^'(r7) follows in a similar way. The proof of (ii) follows from (i) using Note [231 with 
the partition 7V = (0, 1/2, 1). D 

Corollary 4.18. Let Wlwf,W2]{r7) = wf{rY)w2'{rY) — W2{r7)wf'{rY). We have 



sup 

A6[0,l] 



W[wf,w^] ^48-10"^ (91) 

Proof. We substitute A = 1 + 0"^ in ([63]) . and use the polynomials in ([90]) to calculate the 
Wronskian within an accuracy of ±4- 10~^ (obtained by crudely bounding away the effects 
of the errors on the right-hand side of ()90p ). For estimating the resulting polynomials, 
we write o" = ^ + ^z, re-expand and use Note [231 D 



4.1.8. End of the proof of Proposition \4.2\ ' the Wronskian of the actual solutions. 

Lemma 4.19. The Wronskian W[ui,U2]{r7) of the actual solutions satisfies 

\W[ui,U2]{r7) - W[wt,wt]ir7)\ ^ 5 • 10"^ (92) 

Proof. This follows by estimating VF[tii,ti2] via ()9ip . This is straightforward and uses 
Corollarv 14.111 Lemma 14.121 Corollary 14.161 and Lemma 14.171 (ii) to bound |ni 2 — ''^i'2l' 

\u'i2 — C^r^i^2l ^^ ^^^1 ^^ 1^1,2!) 1^2!' 1^1^21 ^^^^ \dr'Wi2\- D 

5. The Operator L_ 

The second main result of this paper is the following one, which establishes the gap 
property for L„ . 

Theorem 5.1. The operator L„ has no eigenvalue or resonance for A in the interval 
(0,1] mL2(M+). 
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As in the case of L+, there are two solutions yi{r; X) and y2{f;X) of the equation ([5]) 



with the properties yi(0;A) = 1 and y2{r;X) 



r e 



l„-rVT^ 



(l + o(l)), r 



oo. Let 



W[yi,y2]{r; A) = yiy'2~y2y'i be the Wronskian of these two special solutions. Theorem 15.1 
is a corollary to the following result. 



Proposition 5.2. One has the lower hound 

sup X-^\W[yi,y2]{rr,X)\ ^ 17/1000 

AG[0,1] 

Therefore yi and y2 are linearly independent for all X G (0, 1] . 
5.1. Proofs. Let 



(93) 



ui(r) = A^i 



yi[r) 



Qjr) 
Q(0) 



; Mr) = J 



y2{r) - -j-Q{r) 



where Ai = A + hi (94) 



(see Lemma l3.10p . Then ui satisfies 

- u" - 2u'/r + (1 - A - Q'^)u = Q/Q{0) 



(95) 



We construct a pair of functions that agree with yi and 2/2 within relatively small errors. 
To this effect, we first define g^2 ™ the following way. Consider the piecewise polynomials 



U,k,l)&S 



'^'; ~92ir) 



kj 



U,k,i)es 



(96) 



where S = {{j,k,l) \ 1 ^ j ^ 3, ^ /c ^ Mj, ^ / ^ 15}, dki-j,eki-j are given in the 
appendix, Xj^ 3 = li2,3, are the characteristic functions of [0, rs), [r^,ri) and [ri,r7] 
respectively, Mj ^ 15 and z depends on r as specified in the top rows of the tables in the 
appendix. To ensure that g'^2 ^^^ C*^ we next let 



?-(r) = 57(r)X47+ ~9-{r) + ~gj{l)-~9^{l-) + {w'^{l) - wAl-)){r - 1) 



Xl4 



+ 



~9-{r)+9j{ri) - ~9j {n-) + 5{g'^{r,) - ~g]{ri-)){r'' - rl)/3\ Xm (97) 
where j = 1,2. We let gj {r) =g'J {r)/r^~'^ and u;j'(r) = {gi {r) - gi {r)\\=Q) / XU 
Lemma 5.3. The following hounds hold: 

|w;j;(r)| ^ 1/100 + 3r/25 V s^ r ^ ry 

r\9i{r)\ ^ TifXoi + '"wXiY, r|52^(r)| ^ %Xm + infXu + ^^47 
\9i\r)\ ^ i§Xo4 + ^47; A92\r)\ ^ f X04 + 3^^X47 



(98) 

(99) 

(100) 



Proof. We check this using Note l4?3l and the partition tt = (O, ri, |,r4,r5, 2,r7). For g^ 
in [0,ri] we first divide ([99]) by r. D 



Define 

and 



i^i 



-w. 



Iwl'lr + (1 - A - Q2)^- _ g/Q(0) 



i?0 = -vj~" - 2w-{'/r + (1 - A - Q2)^- _ Q/Q(0) 
Lemma 5.4. We have \Rq\ ^ 10-^(3x01 + 2^17) and \Ri\ ^ 10~^(6xo4 + 3^47) 



(101) 
(102) 



Note that for L+ we constructed w'^{r) using a difTerent piecewise representation due to the high 
accuracy required. Here, however, it is sufficient to use g^ (r) to define w^(r). 
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Proof. Note that 

l^il ^ \Ro\ + \{Q^ - Q^)w^\ + \Q- QI/Q(o) 

We start with Rq , for which we use the same idea as in Lemma 14.51 (we kept the same 
notations although the functions are different; this should cause no confusion). Due to 
the the monotonicity of Q and the fact that RQ{r)/Q{r)^ is a polynomial we have 

\Ro{ikiz))\ < \Romz))Q\ik{-l))/Qmz)f\ < 3 • 10-5 
using the method in Note 14.31 and the partition given by tt = (O, ^, ^, ^). 



Similarly on [ri , r-j] we use the partition tt 
obtain 



[n 



2 3 4 9 

' 5' 5' 5 



J. 1 3 19 11 „ „ \ „,^j 
' 10 ' ' 2 ' 10 ' 5 ' "' ' / "^^^ 



|i?o(4(^))| < |4(^)/4(-i)Q'(4(-i))/Q(4W)'iio(4(2))| ^ 2- IQ-^ 

The rest of the proof is relatively straightforward, and the details are given in Section [6.41 

D 

The functions g^ , §2 solve a second order equation 

- / + {A{r) - 2/r)g' + {B{r) + A{r)/r)g = (103) 

where 

^^^^ ^ ^:2{r)K"{r)-g{{r)g,"{r) ^^^^ ^^^^ ^ g^ ' {r)^^ " jr) - g^ ' {r)g^ " jr) 
g2ir)g^'{r)-g^{r)g^'{r) g^ ir)g^'ir) -g' {r)g~'{r) 

Lemma 5.5. We have the following bounds on [Ojrj] 

\A{r)\ ^ 2 • 10-4(X[o,^)(0 + ^X[f,rr](^)); \B{r) - 1 + A + Q\r)\ < 3/2500 (104) 

Proof We use the partition TT = (O, ^, |,ri, ^, |, i, |, |, ^, §, 1, i§, |, fi, ^,r6, ry) and 
as usual maximize the numerators and minimize the denominators. 



D 



Note 5.6. We now look for a nearby actual solution ui = w-^ — 6, 5{Q) = 0. In the 
following we write \\f\\ for sup[o^r.7] |/|. 

By dM]) and (fTOTD 5 satisfies -5" - ^6' + (1 - Q^ _ x)s + i^^ = 0, or 

- 6" + {A{r) - 2/r)S' + {B{r) + A{r)/r)5 = /i2(r) (105) 

where /i2(r) = Ri{r) + A{r)6'{r) + B{r) - 1 + A + A{r)/r + Q{rf d{r). 



Lemma 5.7. We have 
\h2{5,r)\ ^10-^ 



| + 2||5'|| + (20 + - 
5 " r 



— + 2||(5'|| + 15| 

10 



Xoi+(5+2||5'||+22|| 



^N,^)+U +""'"+'" 



Xu 
X,n 



Proof. We write 

h2{r)^\Ri{r)\ + \A{r)\\S'{r)\ 



\B{r)-l + X + Q{rf\ 

+ \Air)\/r + |Q(r)2 - Q{rf\\ \5{r)\ (106) 



The estimates now follow from the corresponding bounds satisfied hy Q^ — Q^, Ri,A and 
i?(r) — 1 + A + (5^(r), which are established in (164p . Lemma [5. 41 and Lemma [5.51 On every 
subinterval we use the monotonicity of the exponential and replace it by the value at the 
left endpoint. On the intervals [ri, r4], [r4, 11/5], [11/5, r7] we simply replace 1/r by its 
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left endpoint value. The resulting rational functions are of low degree and are maximized 
explicitly. D 

As in the L_|_ case 5 satisfies 

6 = Hoi6), 6' = H',{S) (107) 

where Hq , Hq are of the form given in (j73p and (I74p with gj' replaced by g~ . 
Lemma 5.8. (i) We have 

No := \Ho{6){r)\ + \Ho{dnr)\/5 ^ 5 • 10"^ + ^^^ + 5 • W-^\\S'\\^ (108) 

Consequently, {Hq,Hq) is a contraction in the ball 

X:={feC\{0,rr))\\\f\\x^6-lO-''} 

where \\f\\x '■= s'^Pr£[o,r7) (1/(^)1 + il/'('")l)- Thus there in an actual solution ui within 
6 • 10~^ in II • ||x of Wi . 

(a) \w'^'{ri) — ui{r-j)\ = |5'(r7)| < 6 • 10~^ (cf. again footnote \3.7\ on p. [T^P . 

Proof. Since 1/W{r) := {g^ {r)g2'{r) — g^'{r)g2{r))~^ = r^exp ( J^' A{s)dsj, we have 

l/|VF(r)| < SlrVsO 

The proof follows from a piecewise analysis of Nq. This is done by straightforward inte- 
gration of the polynomials involved in Nq and then maximization of the rational functions 
(whose numerators are of degree at most 5) multiplying \\6\\, \\5'\\ and of the free term 
in the result of the integration. Taking the derivative of the rational functions and re- 
expanding the numerators at the left endpoint of each interval, we see that all these 
re-expanded polynomials have positive coefficients. Thus the maximum is obtained by 
evaluating Nq at the right endpoint of each interval. We provide the intermediate bounds; 
the largest one, from which the result follows, is on the last interval. 



Nq ^ 10-4|(1/10 + 5\\6\\ + 3/10115' IDxoi + (7/10 + 23||5|| + 3\\6'\\)xu 

+ (4 + 170||5|| + 22||5'||)x[.„n/5) + (5 + 230||5|| + 47||5'||)x[ii/5,h} 

(ii) This follows in the same way as (llOSp . by estimating H'q on [11/5, ry]. D 

Note 5.9. We have ui{0) = w^ (0) = ^5 (A), a quintic polynomial. Substituting A = 
1/2(1 + z) and using Note [231 to estimate this quintic, we get 

|ni (0)1 < 2-10^^ (109) 

Thus Q/Q{0)+Xui is close to but not exactly equal to the function yi defined immediately 
after Theorem 15.11 We will modify it multiplicatively to make up for the discrepancy. 

Lemma 5.10. The function 

n,{r) = n,{r) - u,{0) ^T^^^) ^''°^ 

is the solution of ([95]) with -111(0) = (and thus Q/Q{0) + Xui = yi). 

Proof. Clearly, ui(0) = 0. Using ^, it is straightforward to check that ui solves (|95|) . D 

Lemma 5.11. One has the estimates {w^ {rf)\ < ^ and |u'^'(r7)| < ^. 

Proof. The bound for w^ (rj) follows from Lemma 15.31 We estimate the quintic polyno- 
mial w^'{rj) in A by taking A = 1/2(1 -|- z) using again Note [231 D 
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Proposition 5.12. The solution ui satisfies 

Iw^irj) - ni(r7)| < 7 • 10"^ {w^'irj) - Ui{r7)\ < 7 • 10"^ 

Proof. Note that Lemma 13.61 and the definition of Q imply \Q[rf)/Q[0)\ < jSq and 
|(5'(r7)/Q(0)| < Y§o, and (fT09l) and (fTTOjl imply the (crude) bounds {uiin) - ni(r7)| < 
3-10^^{^ + \ui{r7)\) and \u[{rr) - u[{r7)\ < 3-W~^{-^ + \u[{r7)\). The rest is straight- 
forward from Lemma 15.81 and Lemma 15.111 D 



5.1.1. The region r > r^: the quasi-solution bounded for large r. Let a = \/l — A, and 
ai = (6839/2521)2, ai = 1 + a, a2 = 1 - cr. We consider 

g-{r) = {ar)-'e-^'-{a - aiEi(-2r) + aie2-ViEi(-2air)); w^ (r) = X'^g^ (r) - g^^ir)) 

(111) 
where 

93o{r) = g3 ir)\x=o = e-'' /r + aig{r) (112) 

Since A = 1 — cr^ we have 

2 — 4r 

i?2(r) := -W2" - 2w2'/r + (1 - A - air-^e-^'')w^ - c/30 = "^^ 3 R2 (113) 

where 

R^{r) = ((7a2)~^e-""(2ae(=^+")"Ei(-4r) + (e2^-ae"i'^)Ei(-2r)-e2"iViEi(-2air)) (114) 
Lemma 5.13. We have \R2{r)\ < ^e~'^^'/r^ and \w2 {r)\ < %e~'^'" . 
Proof. Rewriting Ei(— r) as e~^ J-00 ^**/(^ ~ '") ^•^ ^^ obtain 

'''^')-' J_J a^is - Ar){s - 2r){s - 2a,r) 

We use the following inequalities to bound R2{r): For a € (0, 1), r > 0, we have {e"'^^—l) ^ 
(726^; furthermore, s < 0, and thus \s — 2air\ ^ \s — 2r\ and \r/{s — 4r)| < 1/4; we are in 
the range r ^ ry, thus e^*" ^ e^**^, and with m > 0, \s — mr\ ^ \s — mrj]. Therefore, 



1^2 (r) I ^ 







s'^e' 2(2 + e-'^^)se'' 



(s-10)(s-5) 



+ 



4(s-5)^ 



ds = -4e^"Ei(-10) + 2e^Ei(-5) 



-ri , „ , 3a? e- 5 e-2^' 3 e^^r 



- 3Ei(-5)e^^ - V ^ l«2| < ^4^^-^ < :5^^-ir (^5) 

100 r'^ 250 r-^ 

Using (llllj) and (|114p ) we get re^c7iw;^(r) = (e'^^^ — 1)/(T2 — aie^^R2- Therefore, since 
V(T2?' ^ we have e'^^'" — 1 ^ a2re"'^^ we obtain, using the estimate for R2 in (|115p . 

\w, (r)| ^ — < \w, (r)| ^ — ^ ^^ (116) 

and we are done. D 

Since Q satisfies ([3]) and y2 is a solution of ([5]), the definition of U2 (see ([M|) and (fT9|) ) 
implies that U2 satisfies the equation 

- U2" - 2u2'/r + (1 - A - Q'^)u2 = Q/Ai (117) 

Writing U2 = W2 + Sq, we obtain from ()113p and ()117p after regrouping the terms, 

- 6'^ - 26'Jr = -i?2 + (A - 1 + Q2) ^^ + ("gS _ a^^^ ^2" + £ " 53~o =: ^3 (118) 
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Looking for exponentially decreasing solutions, we write (jllSp in the integral form 

/°° tit r°° 
-J s^f{s)ds (119) 



Using (^ and Note SH we get 

-2r 



Q\r)-ai- 



^ 



Qiry - QirY 



+ 



-2r 



Ol- 



Q{r) 



< 



3 e 



-2r 



200 r2 



(120) 



Using Lemma [3.10l to estimate Q in terms of Q and to bound 61^2 ) Definition 12.31 and (jlSp 
to estimate Q, ([10]) to estimate g, ([20]) for A and B, ([94]) for Ai and (fTT2]) we get, for 
r ^ 5/2, 



Q(r)/Ai-r73o(r) 



B5(r) + 62(r) 



^1 



aig{r] 



€ 



b2{r) 



Ai 



+ 



B 



ai g{r) 



< 



1 e 



-3r 



20 r3 
(121) 
In the same way it is shown that for r ^ ry we have Q^ < -^ which implies |1 — A— Q^l ^ 1. 
Therefore denoting ||/|| = sup^^^^ r2e^''|/(r)| we use Lemma [5.131 (|118p . (|120p and 
([m]) to obtain 



|/i3(r)| < (— + 



-2r 



(122) 



Lemma 5.14. Hi is contractive in the ball {/ | ||/|| ^ 1/125}. Furthermore, 1^2(^7) — 
«;^(r7)| < 10"^ and \u'2{r7) - w^^'irj)] < 2 • 10^^ 

Proof. We use the crude estimate |Ti(exp(— 2r)/r^)| < r""^ J^ dt J^ exp{—2s) ds = | 
and obtain from (jl22p 



Iffi(*)l(r)<^(5^ + ^ll*l 



We simply bound |if^(5)| = J^{s/r)'^h^{s) ds using |i7j(exp(— 2s)/s^ 
get 

I^U^o)(.)l^(^ + ^ll^o||)e-/.^ 



(123) 



D 



l^-2„-2r 



r e and 



(124) 



It follows from ^31 and ([IM]) that \So{r7)\ = \Hi{5o){r7)\ < 10-^ |<5^(r7)| = \H[{6o){r7)\ < 
2 • 10"^ 

5.1.2. The Wronskian. The formulas for Wn and Wo' = dj-w^ are 



""^^"'^ = ^r ^'''" " ^) " ^Ei(-10)e^^ 



^^(e'^^.r _ ^)Ei(_5) + ^e^^-Ei(-5ai) (125) 
5A(T 5Acr 



w^'irj) = --|-(2e-^^'^ - 7e-^^ + laie-""'' + bae-'^^'^ - 2aie-''-'n " ^Ei(-10)e'^^ 
25A 25A 

+ 2«iEi(-5) ^^^^_,, ^ 2e-'v- + 5ae-^r^) + ^e'^^-(5a - 2)Ei(-5ai) (126) 
25Acr 25Acr 

It is useful to estimate these in terms of polynomials. 
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Lemma 5.15. (i) With z = 2a — 1 we have 

61 139 97 z2 124 z3 199^4 73 ^5 j^ ^6 ^6 z'^ 24 z^ _ 

560 ~ 588"^"^ 316 22 ~ 40923^ 78624 " 38325"^ 52626 ~ 173 2^ "^ 385 28 ~ ^^ V7, 



23 52 



3 z2 



13 z^ 



19 z4 



22 z5 



21 z^ 



13 z^ 



59 z 



8 



303 587 103 22 31123 235 2^ 257 2^ 292 2^ 242 2^ 1563 2^ 

(a) We also have \w2 {r)\ < 1/2 and {102' {r)\ < 1/5. 
Proof. We consider 



W2 '{rr) 



29 



< 3-10-'' 

< 2-10-4 
(127) 



Ws [(72) 

and 



5 5 



rr . 



'^-^— {e"^'' - a)Ei(-5) + l^e^^"Ei(-5ai) 
5 5 



w^ia) = - ^"f ''' (e"^^^ - a)Ei(-5) + ^e'^^'^Ei(-5ai) 

0(Ti OCTi 



A direct calculation shows that W2 {rj; a) = \ ^Wg {(T2)+cf ^w^ [a] = a^ Jq w^ '(<^2s) ds+ 
L ■w^'{as) ds. The rest of the proof is the same as that of Lemma [4. 171 The calculations 



for ■W2 ' are similar. 



n 



Note 5.16. To estimate the Wronskian W[yi,y2]{r7;X), we express it in terms of ui^2 
using (i94l) . 






Ai 



Q(o) 

+ A (ni(r7)W(r7) - Mr7)uiir7)) (128) 



At the first stage, we look for a nearby quantity solely containing polynomials with rational 
coefficients; we find such an approximation by replacing in (jl28p ui^2 by ^^2' Q by Q, 
and Ai by A (cf. dM]))- Let thus 



^-l^r(^^ .^ ^^1 {r7)Q'(.rY) - w^ [rrYQirr) ^(r?)^^ 'jrr) - Q'{r7)w2 (ry) 
^ '' A Q(0) 

+ A (w;f (r7)u;-'(r7) - u;^(r7)w;-'(r7)) (129) 

Corollary 5.17. The following lower hound holds: 



sup 

Ae[o,i] 

296 



\-^w{\) 



:^ 



55 



(130) 



Proof. By (j20p . we have 1/^4 = gg^ + ^2 where IA2I < 10 . Using the polynomials in 
Lemma I5.15| we can write Vl^ as a polynomial with explicit rational coefficients plus a 
remainder whose absolute value is smaller than 2 • 10~^. Then we re-expand the resulting 
polynomial approximation of M^ at u = 2 using Note 12.51 D 

5.1.3. The Wronskian of the actual solutions. 

Lemma 5.18. The Wronskian of the actual solutions yi,y2, W[yi,y2]{rj; X), satisfies 



\W[yuy2]ir7-A)-W{X)\< 



A 
1000 



VA G [0, 1] 



(131) 
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Proof. The difference between (|128p and (|129p (we omit its long and clearly straightfor- 
ward expression) is bounded simply using the triangle inequality. The terms that need to 
be estimated in the difference are ui^2 — 'Wi2 ^^^ ^i 2~9r'Wi2 ^^ which we use Lemma [5^81 
(ii), Proposition 15.121 and Lemma [5T4| Q/Ai - Q/A and Q/Q(0) - Q/Q{0) which we 
estimate using (pOj) . Lemma [3^ and Lemma [3. 101 I^r2l ^^*^ l'9rW^]~2l which are bounded 
in Lemma 15. Ill and Proposition 15.121 D 

5.1.4. End of the proof of Proposition \5.2l Proposition 15.21 follows from Corollary 15.171 
and (llSip . Therefore W[yi, y2]{rT, A) / for A G (0, 1], implying Proposition 15.21 

6. Appendix 1: Further details of proofs 

6.1. The general solutions of dH) and ([5]). We have shown existence of two solutions 
vq and i^oo of dH (and separately of i^) which belong to L^ near r = and r = oo, 
respectively. Here we address the question of uniqueness. 



Lemma 6.1. Let uq and Uoo be the solutions of (j4]) described in Sections \4-1.3 and \4-l-^ 
Let V be any solution of dH). Then either the Wronskian of v with respect to uq (uoq) is 
zero or else v is not in L^ near zero (infinity, resp.). 

Proof. The functions f/o,oo := fUQ^ooif) satisfy the linear equation 

-U" + {l-X- 3Q^)U = (132) 

The Wronskian of any two solutions of ()132p is constant (since the coefficient of U' is zero). 
If the Wronskian of v with respect to uq (uoo) is not zero, then obviously the Wronskian 
Wq (Woo) of Vi = rv with respect to Uq (Uoo) is a nonzero constant. Now it follows from 
the expression for wf2j Corollarv l4.1H ([86]) . and Lemma [4.151 that Uo{r) — )• 0,C/Q(r) — >• 
ui{0) / as r — ;• 0, and Uoo{r) — )• const, [/^^(r) — ;• {U^ decays exponentially) as 
r — 7- oo. Thus if Wq ^ then either Vi{r) — t- const 7^ or V({r) ~ const/r as r — t- 0, and 
if Woo 7^ then Vi{r) must increase to 00 either exponentially or like const. r as r — ?• 00. 
The conclusion then follows. D 

Lemma 6.2. Let yi andy2 be the solutions of (J5]) described in the beginning of Section\^ 
Let V be any solution of ^. Then either the Wronskian of v with respect to uq (uoo) is 
zero or else v is not in L^ near zero (infinity, resp.). 

Proof. Essentially the same as the proof of Lemma 16.11 using ()94p , Lemma 15.81 ()119p , 
and Lemma 15.141 D 

6.2. Detailed proof of Lemma 13.81 In the sequel, we use the bounds in Lemma 12.41 
(EyI), dMD, dM]), dSHD and Lemma[331 

Note 6.3. We use various bounds to find estimates for the integrals in terms of exponen- 
tials and/or polynomials of low degree, avoiding special functions: (i) We replace co (s), 
s~^ and l/(s -|- 1) by their left endpoint values and ^q:^ by their values at the right end, 

¥^ + ^ ^y ^ ~ 1'^ °^ I-*-/^' ■*■]■ (^^) ^^ ^^^ *^^ partition tt = (i, |, 1) (chosen so that 
the coefficient of the 4th power is small enough not to alter the desired estimate) for 
some quartic polynomials (iii) A function / of the form quadratic plus an exponential is 
estimated by explicitly finding the zeros of the second derivative and inferring the mono- 
tonicity properties of /. (iv) We also note that Hi 2 are different from H[ 2 only by the 
prefactor function. The estimates of H[ 2 are thus obtained easily from those for Hi^2- (v) 
In some integrals where u}{s) would introduce exponentials in the final result, we simply 
bound uj{s) by the value at the left endpoint. (vi) Polynomials of the form x^q{x) with 
q quadratic can be clearly maximized explicitly. 
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6.2.1. r G [0, 1/2). We start with H2 and H2. We break the integral at ^, 1 and | (^ is 
introduced for better bounds; it is not otherwise a special point; in the second integral 
we bound uj{s) by 2e~'*/3). Using Note 16.31 (on [^, 1] we use (v)) to bound the integrands 
we get 



\H2 (S) (r)| ^ f / (1/4 + 3s/20) [pi (11/10 - s) s + Eq (3s/50 + 1/165)] ds 
+ f / (1/4 + 3s/20) [pi (11/10 - s) s + feoe"'' (3s/50 + 1/165)] ds 



+ f / ^-'[P2&-s)s + eo{^, + i,)e-^]ds + f I e-Hps + 2Aeo)^ds 



(133) 



The right-hand side of (jl33p is a quartic polynomial bounded above and below by cubic 
polynomials, obtained by noting that, in the present interval, r^ € [0, 2'^^)- We obtain 

\H2 {6) (r) I < 10"^ (23/25 - 6rV5) + eo/15 (134) 

by explicit extremization of the cubic polynomials. Using Note [631 (iv) we get 



1^2 WMK 



36 
13 



10" 



23 6r' 

25 ~ T 



+ eo/15 



(135) 



To estimate Hi we replace uj{s) = e *(s + 1) ^ by the upper bound 1 so that the integral 
evaluates to a polynomial; we get 

Ift (*) (0 I < (i + I) 5 ^' 13, (., (ii - „) „ + .„ (I + ^) ) ... 

After evaluating the integral, the part without eo is maximized using Note [6]3] (vi) . For 
the derivative we use Note 16.31 (iv): 



\H[{6){r)\<10-'r(2-^^ + '^ 



6.2.2. r G [i 1). Using Note [63] we get 



,,,,,,,,, ^ 7e/ 10 rVl 3 \ / /ll \ /3s 1 \ 2e"i/M 

\H2{6)(r)\ ^ —{ — / - H s] { pi{ s s + eo \ ds 

I 2 w V ; I ^ ^Q I g / \ ^ 20 / r V 10 / V 50 165 / 3 i 



+T /^ '~' {'' (f - ^) ^ + ^° (I + ik) ¥) '^+r '"^ ^^^ + '''°^ ^ '^ 



< 



7e 

To 



10 



-6 



11 

50 



1 



+ 2(l-r))+— eo) (137) 



and 



17/^5) (r) K 18 ( 10-^ ( ii + 2 (1 - r)) + ^eo 



(138) 



32 OVIDIU COSTIN, MIN HUANG, AND WILHELM SCHLAG 

For Hi we distribute e* inside the integral and then use Note 16.31 We get 

For H[ we get 

\H[ (*) (,.) K J3 (s . 10-« + i.„) < (I - >) (s . 10- + l.„) (140) 



6.2.3. r G [?^4,r7). Here we use Note 16.31 except for the re-expansions; other estimates are 

il -10-6 + 10-2(6 _n 

or r^ / / Qo 1 \ ^-1 



explained below. Let ^i = jA • IQ-^ + lO^^^g _ l^r)eo. We get 



/■oo 4e-3* 

+ eW^ e-^(p3 + 24eo)^^ds<Si (141) 

The coefficients of Eq , n = 0, 1 of /i — Bi have the form in Note l6.3l (iii) and are estimated 
as explained there. A nearly identical calculation yields 

\H'2{S){r)\ < |i ^ ^ (142) 

2r r + 1 

In evaluating Hi, we break the interval of integration as follows: [0, g], [^'l] ^^'^ [Ij'^]- 
To further simplify the result, in the integral on [gjl] we first replace s/{s + 1) by its 
(alternating) Taylor series at s = | mod 0{s ^^ 



2' 
2) ■ 



^ I + S(- - 3/4) -^is- 3/4)2 + 256 ^^ _ ^^^^, 



s + 1 7 49 ' ' ' 343 ' ' ' 2401 

On the same interval, we also bound uj{s) ^ 36"^. In the integral on [l,r], we first bound 
s^ /{s + 1) by 5s/7 (which holds for s ^ 5/2) and then apply the bounds in Note [6l3| for 
example, we use that oj{s) ^ 2^~^ ■ After evaluation of the integrals the result is of the 
form described in Note 16. 3|, (iii) and we get 

r\Hi (6) (r) K ^ • lO'^ (3 - r) + ^^o (143) 

Using Note 16.31 (iv) we get 

l^l('^)(OI<'(riO''(3-^) + IM'o) (1^4) 

6.2.4. r ^ rj. The bounds for H2 and H2 are straightforward: We get 
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In the estimate for Hi we proceed as in ^6.2.31 (iii) above, except for the interval [l,rj] 
where instead of the Taylor series of s/{s + 1) we use the simple inequality -^^ ^ M ~ !§• 
After integration, the estimates are elementary and we get 



IH, (6) (r)U- (10- (^- - 36.-J + ^- - -e--J .,J (147) 

6.3. Proof of Corollary 13. 9L We first note that, by definition, Hq ^ 0. Throughout this 
proof, we use the bounds (jl34p ~ (|137p : sometimes we replace them by nearby polynomials 
with simpler coefficients, or majorize them by close enough lower order polynomials, easier 
to maximize. 

6.3.1. r G Ji := [0, i). Here we write 



2' 

1 
2 

-6 



iJo (5) (r)<i3i(r):= 10-5(1 + |r-r2) + Aeo V sC r < - (149) 



Since (r + l)Bi < 35 • 10 in Ji, we have 

I (r + 1) e'^Ho (S) (r) | ^ SSe^/^ . ^q-g < 7 . iq-5 ^^gQ^ 

6.3.2. r e J2 := [i,l). Here 

Ho{6){r)<B2{r):=10-^(—-—r] + (—-—Aeo V-^r<l (151) 

and thus (r + 1) c^Hq ^ maxj2 (r + 1) e'' B2{r) < Eq. Note that the extrema of (r + 1) e^B2{r) 
can be found explicitly. 

6.3.3. r £ J3 := [r^jry). Here we get 

(r + 1) c^Hq (6) (r) < /(r) where /(r) is given by 

lO-^e'' / , /63 21r2\ /91 571 176r2 3, ^2A A , . 

^10-3 +6r--- + - + — r-^ + -(r-l)MeoUeo (152) 



Vio ioyv5 5 5 2 

We artificially added the term |(r — 1)^ to /(r) so that f'{r) is of the form e'^r~'^P^{r) 
with P3 a positive cubic polynomial on J3, as it can be checked by studying its derivative. 
Thus the maximum in (I152p . which is < Eq, is attained at r = ry. 

6.3.4. r £ J = [rj, 00). In this region, (r + 1) e^Ho {6) (r) is bounded by 

^^ . /99 231 33 \ ,^ ^ 2rf 108 126 404 113 \ 

10" 1 1 + 10" e" i 1 

V25 50r 50r2y \^ 25 25r 25r2 lOOr^y 

/ 21 49 7 „2r ( 23 67 431 3 \\ 

"^ ^° V 125 ^ 250r ^ 250r2 "^ ^ V 10 ~ 25^ "^ 50r2 "^ 5^) ) 

. /257 6 107 \ , 2r f 112 26 84 59 \ 

< 10-5 TTT + - + T^T^ +10-^6-2^ _ + + ^^ =: i,+t2 < £0 

\ 50 r 125r7r/ \ 25 or orjr bUr^r J 

(153) 

where the first inequality above follows after replacing eg by 7- 10^^, simple term- by-term 
comparison, and then using the fact that r ^ r7 in this region. For the very last bound 
we first check that ^2 < and then show that (ii — £o)/t2 > — 1 by finding the maximum 
(> —3/10) of this ratio by an elementary computation. 
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6.4. Details of the proof of Lemma l5.4l Combining the result for Rq with Lemma[27 
Lemma 15.31 and Proposition 13. 5| we get 



\R ^^^| < 10-5/3 + 73(1/100 + 3r/25)e-3^- + (7/4)e-7(l + r), rG [0,r2) 

' ^^ ^' ^ \2 + 80{l/W0 + 3r/25)e-'^^''/^/{l + r) + {7/4:)e-''/{l + r), r£[r2,rj) 

(154) 
On the first interval we use the estimate e ■''" ^ 1/(1 + r)-', j = 1,3 to obtain a rational 
function whose numerator is a linear function and whose denominator is (1 + r)^. This 
linear function is estimated in an elementary way, using its derivative. On the second 
interval, since the derivative of (l/100+3r/25)e-^^'''5 is negative, the function 80(1/100+ 
3r/25)e-^^^/5/(l + r) + (7/4)e-''/(l + r) is decreasing. Its upper bound on the interval 
[r2,r4] is 3 and on interval [r4,r7] it is 1. 



7. Appendix: Polynomial quasi- solutions for the soliton 
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23 
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53 
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94 
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600 
19 
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12 


15 
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23 
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J2 
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1106 


3554 


268 


5785 


1239 


5816 


719 


4122 


8351 


21281 
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ON THE SPECTRAL PROPERTIES OF L± IN THREE DIMENSIONS 
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Appendix: Polynomial quasi-solutions for L^ and L_ 



Cki-j for r € [0,r2),z= ^ 



Cki-j for r £ [r2, rg), z = r - 1 



Cfci;j for r e [r5,r7], 
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-1 
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^fc i;;; f"^' "^ ^ [tj, rj], Z = r ■ 
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3 
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efcj.y for r £ [0, n), z ^ Sr 



e-ti-j fo^' '' *= [ri,rA),z 



^kl::, fo'' ^ ^ [?-4,r7),z = r 
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